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ABSTRACT 


In  this  work  a  number  of  semirigorous  and  numerical  techniques 
are  presented  for  analyzing  the  microstrip  transmission  line  corfigurations. 
The  methods  of  quasi-TEM  as  well  as  hybrid-mode  analyses  are  presented 
in  some  detail.  A  discussion  of  the  higher-order  modes  in  such 
transmission  lines  is  also  included, 

In  addition  to  analyzing  the  uniform  line  configurations,  the 
problem  of  evaluating  the  effect  of  losses,  end  loading,  etc,,  are 
examined  and  methods  for  handling  these  problems  are  discussed. 
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I.  INTRODUCTION 


The  microstrip  line  as  it  is  found  today  in  microwave  integrated 
circuit  designs  represents  the  second  generation  of  printed  transmission 
lines.  The  original  version,  called  the  stripline,  was  introduced 
around  1949,  more  than  two  decades  ago!  As  early  as  1955,  a 
special  issue  of  IRE,  published  in  March  1955  as  Transactions  on 
Microwave  Theory  and  Techniques,  was  devoted  entirely  to  the  subject 
of  microwave  strip  circuits.  The  topics  covered  in  this  issue  included 
the  analysis  of  strip  transmission  lines,  computation  of  mutual 
ijnpedance  of  coupled  lines,  radiation  effects,  evaluation  of  discontinuity 
effects,  as  well  as  several  other  related  subjects.  Though  much  interest 
was  expressed  in  this  new  form  of  line  for  a  few  years  immediately  after 
its  introduction,  i ts  use  in  microwave  design  did  not  become  very  popular 
until  quite  rc'****Cj.y  when  new,  low-loss  dielectric  and  ferrite  substrate 
materials  became  available.  The  evolution  of  the  stripline  led  to 
several  different  versions  of  the  microstrip  line,  viz.,  open,  shielded, 
and  boxed  lines,  all  of  which  make  use  of  dielectric  or  ferrite 
materials  as  low-loss  substrate. 

In  this  work  we  will  be  concerned  with  the  analysis  of  several 
of  the  above  configurations  of  the  microstrip  line.  We  will  describe 
a  number  of  semianalytical  and  numerical  techniques  that  have  been 
developed  for  analyzing  these  lines.  It  is  interesting  to  point  out 
that  the  conformal  mapping  technique  used  extensively  in  the  fifties 
t*o  analyze  the  stripline  is  not  conveniently  applied  to  the  microstrip 
line,  since  the  latter  is  an  inhomogeneously  filled  structure.  However, 
at-  will  be  shown  in  Section  II,  a  modified  version  of  the  conformal 
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mapping  technique  can  still  be  used  to  derive  an  approximate  but 
accurate  solution  to  the  microstrip  problem. 

Two  other  important  factors  have  prompted  the  development  of 
new  techniques  for  analyzing  the  microshrip  line.  The  first  of 
these  is  the  advent  of  high-speed  digital  computers  that  have  made 
it  possible  to  apply  numerically  rigorous  techniques  to  the  solution 
of  electromagnetic  and  quasi-static  boundary  value  problems.  The 
second  factor  is  the  increasing  use  of  microstrip  lines  in  the  giga¬ 
hertz  frequency  range  where  the  effect  of  dispersion  in  the  line  is 
no  longer  negligible.  This,  in  turn,  requires  a  full  hybrid-mode 
analysis  of  the  boundary  value  problem.  The  method  of  formulation 
as  well  as  the  solution  of  such  problems  differs  substantially  from 
the  TEM  or  quasi-TEM  approach,  which  are  valid  at  lower  frequencies. 

We  will  illustrate  this  point  more  fully  in  Section  III,  where  we 
will  present  a  hybrid- mode  analysis  of  the  microstrip  line.  Finally, 
the  problems  of  evaluating  the  effects  of  end-loading,  discontinuities, 
and  losses  in  the  microstrip  line,  etc.,  also  require  the  development 
of  techniques  that  have  to  be  tailor-made  for  such  problems.  A 
discussion  of  these  problems  is  also  Included  in  the  following 
sections. 

The  topics  covered  in  this  work  have  been  and  still  appear  to 
remain  very  popular  subjects  for  publications,  as  evidenced  by  the 
prolific  number  of  papers  published  in  the  microwave  literature.  It 
is  physically  impossible  even  to  attempt  to  describe  all  of  the 
contributions  that  have  appeared  in  the  literature  during  the  past 
five  years  on  the  analysis  of  microstrip  lines  and  related  problems. 

We  will  only  present  herein  a  discussion  of  some  of  the  representative 


techniques  and  outline  the  basic  principles  underlying  the 
procedures.  The  reader  interested  in  additional  details  on  a 
particular  method  is  advised  to  consult  the  original  publications 
which  are  referenced  in  this  work. 
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II.  QUASI-TEM  ANALYSIS  FOR  MICROSTRIP  LINE  STRUCTURE 
2.1  Preliminary  Discussion 

Having  presented  a  brief  historical  review,  we  will  now  proceed  to 
discuss  a  number  of  quasi-TEM  techniques  suitable  for  theoretically 
determining  characteristics  of  the  microstrip  line.  As  will  be  evident 
later,  these  analyses  are  valid  only  in  the  low  frequency  range  where 
the  free-space  wavelength  is  much  larger  than  the  strip  width  and  the  thick¬ 
ness  of  the  substrate.  A  more  complete  analysis,  valid  for  arbitrary 
operating  frequencies,  will  be  presented  later. 

Figure  1  shows  some  typical  cross  sections  of  microstrip-type 
transmission  lines.  Before  we  proceed  with  the  details,  however, 
it  will  be  useful  to  explain  why  the  quasi-TEM  analysis  is  restricted 
to  the  low-frequency  region  only.  Though  the  comments  are  applicable 
to  all  of  the  geometries  shown  in  Figure  1,  we  will  illustrate  our 
argument  by  reference  to  the  structure  shown  in  Figure  la.  The  electric 
and  magnetic  field  components  of  the  characteristic  solution  in  this 
structure  may  be  expressed  in  terms  of  a  scalar  potential  <j>  which  is 
required  to  satisfy 

+  (k^  -  3^)<£  -  0  in  the  air,  (la) 

2  2  2 

V  *  +  (e^k  —  B  )4>  —  0  in  the  substrate,  (lb) 

where  k  =  2n/A  is  the  free-space  wavenumber,  and  3  is  the  propagation 
constant.  In  the  limit  of  =  1,  the  lowest-order  solution  for  (1) 
is  TEM  with  3  =  k.  For  this  limiting  case,  the  potential  <{>  satisfies 
the  Laplace  equation  in  the  cross  section  and  does  not  generate  longitudinal 


4 


electric  or  magnetic  field  components.  For  the  practical  case  of 

^  1,  Equation  (lb)  may  be  regarded  as  a  perturbation  of  the  limiting 
case  =  1.  The  effect  of  this  perturbation  is  small  when  the 
operating  frequency,  i.  e.,  k  is  small.  The  quasi-TEM  approximation 
may  then  be  regarded  as  a  zero-order  solution,  to  the  exact  equations 
shown  in  (1).  Though  much  simplification  in  the  analysis  results  with 
the  use  of  the  quasi-TEM  approximation,  it  should  be  kept  in  mind  that 
the  results  obtained  via  this  method  are-  hot  accurate  for  larger  values 
of  k. 

Under  the  TEM  approximation,  the  microstrip  line  is  adequately 
described  in  terms  of  two  basic  parameters,  viz.,  the  characteristic 
impedance  Z  and  the  propagation  constant  y.  It  is  well  known  that 
for  a  low-loss  transmission  line 

2  ^  +  2  (jwL  "  jmC ^ 

Y  =  ot  +  jg  (3a) 

a  =  ~  (£  +  §) ,  6  =  (3b) 

where  R,  G,  L,  C  are  the  resistance,  conductance,  inductance,  and 
capacitance  per  unit  length,  respectively,  of  the  infinitely  long 
transmission  line;  w  is  the  angular  frequency;  and  a  and  8  are  the 
attenuation  and  phase  constants,  respectively.  For  a  lossless  line 
R  =  G  =  0,  and  hence,  a  =  0,and  Z  =  /L/C  is  a  real  quantity.  The 
values  of  8  and  Z  can  also  be  expressed  in  terms  of  the  phase  velocity 
v  =  l/v'Ec  as 


8  =  rn/v,  Z  =  l/(vC) . 


(4) 
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It  should  be  mentioned  that  dispersion  effects  are  neglected  in  the  TEM 
approximation;  hence,  the  phase  velocity  is  equal  to  the  group  velocity. 

Next  we  will  show  that  for  the  lossless  TEM  system  the  boundary 
value  problem  associated  with  the  determination  of  E  (and  Z)  reduces 
to  that  of  calculating  the  capacitance  per  unit  length  C.  Consider  the 
two  configurations  of  the  microstrip  line  shown  in  Figure  la.  Note 
that  one  of  these  has  a  dielectric  substrate  with  relative  permittivity 
e^  while  the  companion  structure  is  a  homogeneous  TEM  line  obtained  by 
removing  the  substrate,  i.  e.,  by  letting  =  1.  It  is  easy  to  show 
that  for  the  latter  case 

6o  =  w/c,  Zq  =  l/(cCo)  (5) 

where  c  is  the  velocity  of  light  in  free  space.  The  subscript  0  in 
(5)  serves  to  distinguish  the  line  from  the  inhomogeneous  microstrip 
line  in  which  ^  1.  If  the  substrate  is  nonmagnetic,  the  value  of 
L,  the  inductance  per  unit  length,  is  assumed  to  be  identical  for  both 
lines.  Hence 


where  Zq  and  BQ  are  given  by  (5) .  Equation  (6)  implies  that  Z  and  6 
are  obtainable  from  the  capacitance  per  unit  length,  C  and  Cq,  of  the 
microstrip  and  unloaded  lines. 

Next  we  will  describe  a  number  of  quasi-TEM  techniques  for 
evaluating  the  capacitance  of  microstrip  lines  of  the  type  shown  in 
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Figure  1.  The  discussion  below  will  include  material  from  a  number 
of  papers  that  have  appeared  in  the  literature  since  1965.  As- mentioned 
earlier,  a  more  rigorous  analysis  based  upon  the  hybrid-mode  approach 
will  be  presented  in  a  later  section. 

2.2  Modified  Conformal  Mapping  Method 

This  method  was  introduced  by  Wheeler  (1965)  for  computing  the 
capacitance  per  unit  length  of  microstrip  lines  of  the  type  shown 
in  Figure  la.  The  concept  of  the  so-called  "filling  factor,"  also 
introduced  by  Wheeler,  is  considered  to  be  a  rather  convenient  and 
useful  way  of  describing  the  parameters  of  microstrip  structures. 

The  right  half  of  the  cross  section  of  the  microstrip  line  is 
shown  in  Figure  2  where  a  magnetic  wall  is  placed  along  the  axis  of 
symmetry,  the  y-axis.  The  strip  thickness  is  assumed  to  be  infinitesimally 
small. 

The  first  step  in  Wheeler's  technique  entails  the  application  of 
conformal  mapping  to  the  geometry  under  consideration  with  the  result 
that  a  simpler,  parallel-plate  geometry  is  obtained  in  the  new  domain. 

For  the  wide  strip,  the  mapping  function  is  chosen  to  be 

z  =  jn  +  d  tanh  —  z'  -  z'  (7) 

with  d  =  g'  for  g'/2  >>  1  where  g'  is  the  effective  width  of  the 
parallel  plate  in  the  mapped  plane  shown  in  Figure  3a. 

Note  that  the  mapping  from  z  to  z'  transforms  the  original 
microstrip  line  into  a  parallel-plate  structure  bounded  by  two 
vertical  magnetic  sidewalls  (x'  =  0  and  g').  Also,  the  planar  dielectric- 
air  interface Q)  - Q)  in  Figure  2  maps  into  a  curved  surface@  -Q)  in 


8 


/, 


Figure  3a.  Although  the  geometry  in  the  transformed  plane  is  seemingly 
simpler  to  analyze,  it  should  be  -realized  that  the  space  between  the 
parallel  plates  is  now  inhomogeneously  iiiled  with  a  dielectric  which 
has  a  curved  boundary.  Since  an  exact  solution  is  difficult  to  obtain 
for  the  partially  filled  capacitor  problem,  in  Figure  3a,.  it  is  necessary 
at  this  stage  to  introduce  some  approximations.  In  Wheeler's  approach, 
the  approximation  is  introduced  in  the  following  manner.  Let  the 
area  enclosed  by  the  curve  Q)  -rQ),  and’  the  lines  (3)  -  and  (3^  ~(7)  in 
Figure  3a  be  ns'.  Let  this  area  be- approximated  by  two  rectangles  with 
areas  its"  and  ir(s*  -  s")  as  shown  in  Figure  3b.  The  area  us"  effectively 
adds  s"  to  the  width  of  the  completely  filled  region  on  the  right 
(a1  <  x'  <  g')  and  hence  is  termed  the  "parallel"  portion.  The  other 
area  7t(s'  -  s")  is  effectively  in  series  with  the  free-space  region 
irso.  Combining  these  two  effects  we  can  approximately  replace  the 
original,  partially  filled  rectangular  region  bounded  byCjQO©  in 
Figure  3a  by  an  equivalent  geometry (JX^Q)  ©>  shown  in  Figure  3c,  in 
which  the  dielectric  filling  has  the  width  s.  It  is  obvious  that  the 
capacitance  calculation  is  straightforward  in  the  equivalent  geometry 
due  to  the  planar  nature  of  the  dielectric  filling. 

The  effective  width  s  may  be  expressed  as 


o'  -  q" 


s  =  s"  + 


(9) 


The  effective  filling  factor  can  be  defined  by 


„  _  S1  -  a'  +  s 
q  "  g' 


(9) 


The  effective  dielectric  constant  e  ^  may  now  be  expressed  in  terms 
of  the  filling  factor  q  via  the  relation 
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(10) 


e 


eff 


(1  -  q)  +  qer. 


Finally,  the  capacitance  per  unit  length  is  given  by  C  =  ee£fCo  where  Cq 
is  the  capacitance  per  unit  length  of  the  unloaded  TEM  line.  It  is 
evident  that  the  problem  is  now  reduced  to  that  of  deriving  approximate 
values  of  the  quantities  s',  s",  or  q. 

An  approximate  expression  for  wide  strips  has  been  provided  by 
Wheeler  that  reads 


s'  =  0.732  [a'  -  cosh-1(0.358  cosh  a'  +  0.953)]  (11a) 

s"  =  .0.386  -  l/2(d  -  1).  (lib) 

The  effective  dielectric  constant  e  can  be  derived  by  the  use  of 

eff 

(11)  in  (8),  (9),  and  (10).  Also,  he  has  shown  that  in  the  case  of 
narrow  strips,  the  effective  dielectric  constant  is  given  by 

e  +  1  g' (e  -  1)  . 

eeff  =  ~  2 . + - L_2 -  [Zn  77/2  +~r  Zn  */*]  *  (12> 

2ir  r 


2.3  Finite-Difference  Method  (Relaxation  Method) 

This  method  is  strictly  a  numerical  technique  for  determining  the 
cross-sectional  field  distribution;  this  knowledge  allows  one  to  compute 
the  desired  characteristic  parameters  of  the  transmission  line.  The 
essential  step  in  this  method  is  to  quantize  the  unknown  field  distribution 
at  discrete  intersections  of  coordinate  grids  called  the  net-points  or 
mesh-points.  The  relaxation  method  is  a  numerical  algorithm  for  solving 
the  unknowns  at  a  large  grid  of  the  net-points  by  applying  the  method 
of  successive  approximations. 

Figure  4a  shows  a  portion  of  the  distribution  of  such  net-points. 

Let  the  potential  at  the  net-point  A  be  New  let  the  potential 
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<p(x,  y)  be  expanded  in  Taylor’s  series  about  the  point  A  to  find  the 
potentials  at  the  adjacent  net-pointy  B,  C,  D,  and  E,  The  expressions 
r.aad 


<j)„  =  $,  +  h 


m  +4-  h2 1 

(i5&J 

lax/.  ^2! 

A 

ax2/ 

m  h2 

' 

>  3y  'A  21 

ay2/ 

(li|  +4  h2 

(A 

l3xlA  21  n  1 

ax2, 

mA  h2i 

'  a/-<j) 

.  2 

A 

ay 

(13) 


+  0(h) 
A 


+  0(ti  ) 


where  0(h4)  implies  the  terms  in  h^  and  higher.  Adding  the  above 


expressions  for  the  potentials,  we  obtain  the  equation 


+  *c  +  *d  +  *E  =  4*a  +  h 


+  0(h  ) . 


(14) 


Utilizing  the  fact  that  the  potential  $  satisfies  Laplace’s  equation, 

2 

we  note  that  the  coefficient  of  h  in  the  right-hand  side  of  (14)  is  zero 

4 

Further,  if  we  neglect  the  terms  of  0(h  ),  Equation  (14)  reduces  to 

*A  =  T  (*B  + +  (15) 

The  numerical  solution  to  this  equation  can  be  obtained  by 
systematically  assuming  the  net-point  potential  at  each  of  the  points 
where  the  potential  is  unknown,  numerically  testing  the  residual  error 
in  this  assumption,  and  using  the  residual  error  in  the  following 
way  to  update  the  assumed  potential.  The  pertinent  equations  for 
applying  the  relaxation  method  are 
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,(n)  _  (n) 

A  "  *A 


(n+1)  _  .(n)  (n) 

'A  *  ‘  “rA 


(16b) 


where  the  quantities  with  superscript  n  are  to  be  associated  with  the 

values  obtained  after  the  n^  iteration  and  is  the  so-called 

A 

residual  at  the  net-point  A.  Equation  (16b)  represents  the  updated 

til 

estimate  of  the  net-point  potential  for  the  (n+1)  iteration  and  the 
constant  a  is  called  the  accelerating  factor.  The  iteration  equations 
‘  in  (16)  are  repeatedly  used  until  the  values  of  the  residual  become 
smaller  than  some  prescribed  value  for  all  the  net  points. 

The  iteration  method  is  convergent  as  long  as  0  <  a  <  2.  The 
range  0  <  a  <  1  is  referred  to  asunderrelaxation  whereas  1  <  a  <  2 
is  called  the  range  of  overrelaxation  for  the  acceleration  factor. 

The  convergence  of  the  procedure  is  most  rapid  for  a  somewhere  between 
1  and  2  (Green,  1965)  and  the  iterative  procedure  with  a  in  this  range 
is  referred  to  in  the  literature  as  the  method  of  successive  over¬ 
relaxation  (SOR) . 

The  expression  for  the  potential  requires  a  modification  at  the 
interface  between  air  and  the  dielectric,  substrate  (see  Figure  4b) . 

The  appropriate  equation  to  use  instead  of  (15)  (Green,  1965)  is 

♦a-T«b  +  V  +  2(TTl<er*c  +  V-  (17) 

Finite  difference  equations  at  the  various  boundary  points  have  been 
tabulated  by  Green  (1965)  and  some  general  discussion  can  be  found 
in  a  paper  by  Wexler  (1969).  Interested  readers  are  encouraged  to 
consult  these  publications  for  further  details. 
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It  is  worthwhile  to  mention  that  the  method  of  solution  of 
Equation  (15)  is  by  no  means  the  only  one  available  to  us.  In  fact 
one  can  derive  a  matrix  equation  of  the  form 


by  applying  either  (15)  or  its  variants  [e.g.  (14)  ]:  to  all  of  the 
net  points  in  the  cross  section.  The  matrix  equation  (18)  can  now 
be  inverted  to  obtain  the  solution.  However,  the  disadvantage  of  this 
approach  is  that  a  large  coefficient  matrir  is  required  to  obtain  a 
reasonably  accurate  solution.  Furthermore,  since  a  large  number  of 
the  elements  of  the  coefficient  matrix  are  zero,  the  matrix  is  very 
sparse;  the  direct  inversion  of  Equation  (18)  is  numerically  inferior 
to  the  SOR  technique. 

Once  the  potential  at  each  net-point  is  calculated,  it  is  relatively 
straightforward  to  calculate  the  capacitance  of  the  line. 

One  uses  the  formula 


P  dZ  er  ||  dZ  (19) 

L  L 

for  calculating  the  charge  per  unit  length  where  L  is  the  contour  sur¬ 
rounding  the  center  strip,  n  is  the  outward  normal,  and  is  unity 
if  the  net-point  is  in  the  air  and  equal  to  the  relative  permittivity  of  the 
dielectric  otherwise  .  Equation  (19)  may  bewritten  in  a  discretized  forr 
as 
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(20a) 


o  -  *,<■  i  p. 

S  1=-1  1 


(”) 

WP; 


h  - 


(20b) 


In  this  case,  the  prime  on  the  summation  in  (20a)  indicates  that  the  first  and 

the  last  terms  in  the  summation  are  to  be  halved  in  order  to  approximate 

the  line  integral  by  the  trapezoidal  rule.  The  summation  over  s 

implies  that  the  contour  is  subdivided  into  s  straight  sections.  The 

point  P  is  taken  along  the  dotted  line  (contour)  shown  in  Figure  5 

(Green,  1965).  The  capacitance  per  unit  length  is  readily  obtained 

from  the  value  of  Q  by  using 


where  V  is  the  potential  difference  between  the  strip  and  the  ground 
plane.  To  evaluate  the  characteristic  parameters,  say  the  characteristic 
impedance  Z  and  the  phase  constant  g,  we  need  to  calculate  Cq  for  the 
case  where  the  substrate  is  removed  Equation  (6)  may  then  be  used 
to  obtain  Z  and  6, 

Figure  6  shows  a  few  examples  of  calculated  curves  (Stinehelfer , 
1968)  of  Z  for  the  microstrip  line  enclosed  in  a  shield  (see  Figure 
lb).  For  comparison  purposes,  Wheeler ’ s  results  for  a  flat  (open) 
microstrip  line  are  also  shown  in  Figure  6, 


2.4  Variational  Method  in  the  Fourier  Transform  Domain 


This  method,  which  uses  a  combination  of  variational  technique 
and  the  Fourier  transform  approach,  was  first  introduced  by  Yamashita 
and  Mittra  in  1968  for  calculating  the  characteristic  parameters  of 
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SURFACE  OF-  INTEGRATION 


Figure  5. 


Surface  of  integration  for  calculating  the  total  charge  on  the 
strip  conductor. 


microstrip  lines  and  other  related  structures  (Yamashita  and  Mittra, 
1968;  Yamashita,  1968).  It  has  been  shown  by  the  above  authors  that 
the  method  is  numerically  very  efficient. 

For  the  sake  of  illustration  we  consider  the  shielded  double¬ 


layer  microstrip  line  shown  in  Figure  7.  The  conventional  microstrip 

*  A  A 

line  may  be  obtained  from  this  structure  by  letting  = 

and  d  where  e^,  e,^  and  are  relative  dielectric  constants. 

Once  again  we  will  compute  the  capacitance  per  unit  length  along 
the  uniform  microstrip  line.  To  this  end,  we  first  write  Poisson’s 
equation  for  the  potential  distribution  <j>(x,  y)  in  the  cross-sectional 
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GROUND  PLANE 


Figure  7.  Shielded  double-layer  microstrip  line. 
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where  <K3,  y)  is  the  Fourier  transform  of  <p(x,  y)  and  p(6)  is  the 

transform  of  p(x). 

The  general  solution  for  y^h  +  s+  p  is  a  combination  of 
exp(6y)  and  exp(-6y).  For  the  limiting  case  of  d  ->  ®,  only  the 
exp (-By)  type  of  solution  is  retained  in  the  unbounded  region. 

Returning  to  the  case  of  finite  d,  the  boundary  and  interface  conditions 
for  the  transformed  potential  are  given  by 


?(3,  o)  =  o 

(24a) 

£(3,  h+o)  =  £(8,  h-0) 

(24b) 

ej  h+0)  =  e*-^?(8,  h-0) 

(24c) 

4>(3,  h+s+0)  =  <j>(3,  h+s-0) 

(24d) 

C3dy  <^»  h+s+°)  =  e2  dy  h+s'''5> 

(24e) 

(j>(8,  h+s+p+0)  =  <f(8,  h+s+p-0) 

(24f) 

*  d  <v, . .  *  H  'v  . 

e3  £(e’  k+s+p+0)  =  £*  ~  ?(&,  h+s+p-0)  -  j-  p(B) 

o 


(24g) 

d>  (8 ,  h+s+d)  =  0  .  (24h) 

After  incorporating  these  boundary  conditions  into  the  general 

representation  of  the  solutions ,  we  obtain  a  set  of  linear  inhomogeneous 

equations  for  the  unknown  coefficients  of  potential  functions.  Letting 

p  ->•  0  the  solution  for  the  potential  on  the  strip  is  found  to  be 

?(8,  h+s)  =  —  p  (8)  g(B)  (25) 

o 

where  eq  is  the  free-space  permittivity  and 
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is  the  total  charge  on  the  strip  per  unit  length.  Using  Parseval's 
formula,  Equation  (28)  can  be  converted  into  the  transform  domain 

P  (8s  £(8,  h+s)d8.  (30) 

The  above  formula  is  numerically  more  efficient  to  handle  than  the 
corresponding  space  domain  version  (28),  because  <f>  in  (30)  is  just 
a  product  of  two  functions  while  $  in  (28)  is  a  convolution  integral 
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represented  by  (27).  This  is,  in  fact,  the  major  advantage  to  be 
gained  by  the  use  of  the  transform  technique. 

So  far,  we  have  considered  only  an  infinitely  thin  strip. 

However,  it  is  possible  to  approximately  extend  the  method  to  the 
case  of  finite  t  providing  t  is  not  too  large.  This  is  done  by  replacing 
<J>(S,  h+s)  in  (30)  by  the  average  of  ^(S,  h+s)and  $((3,  h+s+t)  where 
$(S,  h+s+t)  is  the  potential  distribution  at  y  =  h  +  s  +  t.  The 
expression  for  the  transform  of  this  potential  is 

*«•  h+s+c>  ■  * <6’  h+s)  •  (31) 

The  line  capacitance  may  now  be  expressed  as 

00 

£  ■  —K  I  P(8)  »(B,  h+s)  iS(B)  as  (32) 

C  2"Q  -L 


where 


h(6)  =  j  {1  + 


sinhf Is 

(d  -  t)]. 

sinh( 

6  d) 

(33) 


The  final  form  of  the  line  capacitance  including  the  strip  thickness 
is 

00 

i  =  — l —  f  [p(3)]2  gC0)  h(3)  dS.  (34) 

7iO  e  J 
o  o 

The  unloaded  line  capacitance  is  obtained  by  letting 

^  ^  ^ 

=  ^2  =  e3  =  1  in  (34) . 

Although  p(8)  is  still  unknown,  the  variational  nature  of  (34) 
introduces  only  a  second-order  error  when  an  approximate  form  of 
p(B)  is  substituted  into  it.  Since  the  stationary  expression  (34) 
gives  an  upper  bound  of  1/C,  the  calculated  value  of  C  is  always 
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smaller  than  the  correct  value.  Hence  the  choice  of  p(8)  that 
maximizes  the  value  of  C  clearly  yields  the  result  that  is  closest 
to  the  exact  value  for  the  capacitance.  Several  trial  functions  for 
this  problem  have  been  tested  by  Yamashita  and  Mittra  (1968)  and 
Yamashita  (1968), and  numerical  results  have  been  reported  for  two 


different  trial  functions.  These  are 


(m, 

(1)  p  (x)  =  < 

1° 


otherwise 


(35a) 


p(8)  2  sin(6  w/2)  sin(8  w/4)  2 

Q  8  w/2  0  w/4 

/I  +  ] 2x/w| 3,  <  x  <  | 

(2)  p(x)  =/ 

1  0  otherwise 


(35b) 

(36a) 


Figures  8-12  present  the  numerical  results  for  the  microstrip 
^  ^ 

lines  where  =  c  ,  =  1  and  d  •+  ».  The  trial  function  in 

(35b)  has  been  used  for  computing  these  data.  Figure  9  presents  a 
comparison  of  the  results  obtained  by  this  process  with  those  of 
Wheeler.  In  addition,  a  comparison  with  the  experimental  results 
obtained  by  Arditi  (1955),  Dukes  (1956),  and  Hyltin  (1965)  is 
exhibited  in  Figures  10  and  12.  Numerical  results  for  the  shielded 
microstrip  line  were  calculated  using  the  expression  in  Equation  (36b) 
for  the  trial  charge  distribution.  These  data  are  presented  in 
Figures  13-16.  Typical  computation  time  was  less  than  10  seconds 
on  an  IBM  7094  computer. 
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Figure  9.  Calculated  result  for  characteristic 
with  the  results  by  Wheeler  (1965). 
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Figure  15.  The  calculated  guide  wavelength,  e*  =  9.9  (Sapphire);  c*  =  1; 
s  =  0;  t  =  0.  J 


Figure  16.  The  calculated  guide  wavelength, 
s  =  0;  t  =  0.02h. 


e*  =  9.9  (Sapphire);  e*  =  1; 


Before  terminating  this  discussion,  it  will  be  useful  to  present 
a  systematic  procedure  for  improving  the  approximate  solution  for  the 
charge  distribution.  This  is  accomplished  via  the  application  of 
Galerkin’s  method  or  Rayleigh-Ritz  procedure  in  the  transform  domain 
(Ward,  Mittra,  and  Itoh,  1971).  The  procedure  is  outlined  in  the  following. 

Let  us  first  rewrite  Equation  (25)  as  follows 


4>±(C,  h+s)  +4>0<6,  h+s)  =  ~  p(0)  g  (B) 


where 


w 

^(6,  h+s)  =  J 


4>(x,  h+s)  e^x  dx 


(38a) 


<J>o(B,  h+s)  = 


<l>(x,  h+s)  e^X  dx  + 


<j>(x,  h+s)e^X  dx. 


(33b) 

It  should  be  noted  that  is  known  since  4>(x,  h+s)  is  given  for 

|x|  <  w/2  on  the  strip.  However  <fo,  the  transform  of  |x|  s  w/2,  is  not  known 

since  (x ,  h+s)  is  as  yet  unknown  for  |x|  >  w/2  (outside  the  strip). 

The  next  step  is  to  apply  Galerkin's  method  to  Equation  (37).  To  this 

end  p(B)  is  expanded  in  terms  of  known  basis  functions  p  ($)  with  the 

n 

following  unknown  weight  coefficients 


P  (8)  =  l  c  p  (6) . 
,  n  n 
n=l 


We  choose  the  basis  functions  p  (8)  such  that  their  inverse  transforms 

n 

have  finite  support  in  the  space  domain,  i.  e., 


00 

pn(x)  ‘  27  /  K 


(6)  e"j6x  d0  =  0, 
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| x |  >  w/2 . 

(40) 


Next  Equation  (39)  is  substituted  into  (37)  and  the  inner  product  is 
taken  of  both  sides  of  the  resultant  equation  with  Pn(3)>  n  =  1,  2,  « 
This  leads  to  the  following  matrix  equation 


b  =  y  K  c 
m  L.  mn  n 
n=l 


m  =  1,  2, 


where 


»  w// 

b  «  f  p(B)  £,(B,  h+s)dg  =  ~  f  p  (x)4>  (x,  h+s)dx 
m  J  mi  Zu  I  mi 

^  J  tr\ 


CO 

-f  / 

o  ^ 


P  (6)  g(B)  P  (g)dg. 
m  n 


(42a) 


(42b) 


It  should  be  mentioned  that  (42a)  has  been  derived  via  the  use  of 
Parseval's  formula.  An  application  of  this  formula  helps  to  eliminate 
unknown  <f>Q  which  is  no  longer  present  in  (42a) .  The  integral  of  the 

*\i 

product  <J>q  and  does  not  contribute  to  (42a)  because  the  inverse 
transforms  of  these  functions  are  nonzero  in  complementing  regions 
only,  and  hence  their  product  is  identically  zero. 

The  next  step  in  the  solution  is  to  invert  the  matrix  Equation  (41) 


and  solve  for  the  coefficients  cn«  The  line  capacitance  can  then  be 
obtained  in  terms  of  the  c^  by  using  the  formula 


m  w/  i 

C  =  \  c  f  p 

V  L.  n  /  n 

t  n=l 


(x)dx 


where  v  is  the  potential  difference  between  the  strip  and  the  ground 
plane. 

It  is  well  known  that  the  expression  in  (43)  gives  a  stationary 
value  of  C  and  that  the  accuracy  of  the  result  may  be  improved  by 


increasing  N.  Experience  has  shown  that  the  increase  in  the  accuracy 
is  usually  rather  slow  after  initial  improvement  when  N  is  increased 
from  a  very  small  value. 

2.5  Integral  Equation  Method 

So  far  we  have  discussed  three  different  techniques  for  determining 
the  characteristics  of  the  microstrip  line.  Two  of  these,  viz.,  the 
modified  conformal  mapping  and  the  variational  methods  find  their 
principal  use  in  the  cases  where  the  center  strip  is  thin.  However, 
for  some  practical  microwave  integrated  circuit  structures,  the  effect 
of  finite  thickness  of  the  center  strip  cannot  be  neglected. 

Since  the  relaxation  method  is  principally  a  numerical  technique, 
it  can  be  applied  to  the  case  of  the  arbitrarily  thick  center  strip. 
However,  as  mentioned  earlier,  the  drawback  of  this  method  is  that 
its  numerical  convergence  is  rather  slow. 

In  this  section,  an  integral  equation  approach  will  be  presented 
that  can  be  used  to  solve  the  thick-strip  problem  in  a  numerically 
efficient  manner  (Yamashita  and  Atsuki,  1970),  Consider  the  cross 
section  of  the  transmission  line  shown  in  Figure  17.  The  potential 
satisfies  the  Poisson’s  equation. 

?j>(x,  y)=  -  “P(x,  y) .  (44) 

Define  a  Green's  function  G(x,  y , ;  xq,  y  )  via 

V^G(x,  y;  xq,  yQ)  =  -  <S(x  -  xq)  6 (y  -  yQ)  (45) 

where  G  satisfies  the  same  boundary  and  continuity  conditions  except 
for  the  source  condition.  The  Green's  function  is  the  potential  at 
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the  point  (x,  y)  due  to  a  unit  charge  located  at  y  ) .  Applying 
the  superposition  principle,  the  potential  function  4>(x,  y)  in  (44) 
may  be  expressed  as 


<Kx,  y)  =  /  G(x,  y;  xq,  yo)p(xQ,  yo)d£Q  (46) 

where  the  integral  is  defined  over  the  conductor  surface. 

To  find  the  Green's  function  we  proceed  as  follows:  First,  we 
expand  the  Green's  function  in  a  Fourier  series  in  the  x-coordinatet 
Next,  the  partial  differential  equation  (45)  is  reduced  to  a  set  of 
ordinary  differential  equations  of  the  variable  y„  The  solution  of 
each  ordinary  differential  equation  in  each  region  of  dielectric 
materials  is  a  linear  combination  of  hyperbolic  functions-  When  the 
boundary  and  interface  conditions  are  applied  to  each  of  these  solutions, 
the  amplitude  coefficients  of  these  functions  generate  a  set  of  linear 
inhomogeneous  equations  which  are  subsequently  solved  for  the  unknown 
amplitude  coefficients*  Substituting  these  coefficients  in  the  equation 
for  the  Green's  function  completes  the  derivation  of  the  desired 
expression  for  the  Green's  function.  When  the  source  located  in  the 
range  h^  +  <  yQ  <  b  and  the  boundary  condition  at  the  outer  conductor 

is  G  =  0,  the  expression  for  the  Green's  function  is  given  by  (Yamashita 


and  Atsuki,  1970,  1971) 


y;xo, 

: 

1! 

o 

2  n 

n=x 
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o  S 
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2 
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hl  +  h2  1  -yo-  b 


(47) 


where 


*  *  nirh  nirh9 

A  (y)  =  e. e0  cosh  (  - -  )  cosh  (  -  )  sinh  (  — ^) 

n  1  2  a  a  a 

*  *  n7rhi  n7lh9 

+  e0e0  sinh  (  -  )  cosh  ( - )  cosh  (  — *-) 

2  3  a  a  a 

*  *  n7llli  nilh? 

+  e-e,  cosh  (  -  )  sinh  ( - )  cosh  (  — *-) 

j>  x  cl  Si  2 

Ao  nwh.  nuh, 

+  e.  sinh  (  - -  )  sinh  (  - -  )  sinh  (  )  . 

2  a  a  a 


The  representations  of  the  Green’s  functions  for  the  first-  and  second¬ 
er 

dielectric  regions  have  been  omitted  here  as  they  are  not  needed  for 


subsequent  derivation. 

By  following  the  conventional  methods  for  deriving  the  integral 


equation  for  the  charge  distribution,  one  arrives  at  the  desired  equation 


<f>(x,  y)  =  V  »  f  G(x,  y;  xq,  yo)p(xQ,  y0)dfco. 


The  equation  can  now  be  transformed  into  a  matrix  form  by  the  conventional 
discretization  procedure.  The  results  obtained  by  solving  this  matrix 
equation  are  presented  in  Figure  18. 

Typical  computation  time  of  the  characteristic  impedance  per  one 
structure  is  about  30  seconds  when  the  strip  conductor  was  divided  into 
30  sections  (Yamashita  and  Atsuki,  1971). 


2.6  Generalized  Wiener-Hopf  Techniques 

In  contrast  to  some  other  approaches,  the  method  to  be  described 

* 

However  the  interested  reader  may  find  the  pertinent  expressions  in 
Yamashita  and  Atsuki  (1970). 
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Figure  18. 


Characteristic  impedance  of  thick-strip  line  of  Figure  17 

r*  =  rrt  =7.  —  n  it.  t  .  _  ° 


C1  “  e3  =  1;  e2  =  9,35;  hi  =  h3  "  0.4b;  h2  =  0.2b;~w  =  0. 
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below  is  useful  for  deriving  accurate  numerical  results  for  both  the 
charge  and  the  potential  distributions  for  the  microstrip  line  shown 
in  Figure  7«  The  technique  was  presented  in  a  recent  paper  by  Mittra 
and  Itoh  (1970).  The  starting  point  is  the  equation  for  the  charge 
distribution  p  in  the  transformed  domain.  This  equation,  which  was 
derived  earlier  in  Section  2.6,  is  reproduced  here  for  reference 
purposes.  The  transformed  equation  reads 

j-  P (8)8(0)  =  £(8,  h+s).  (25) 

o 

A  semirigorous  rather  than  variational  method  for  solving  this 
equation  will  now  be  discussed.  Note  first  of  all  that  the  potential 
of  the  center  strip  is  one  volt.  Thus,  we  can  write  the  potential 
function  at  y  =  h  +  s,  i.  e.,  in  the  plane  of  the  strip,  in  the  form 

4»(x,  h+s)  *  h(x  -  |)  U(x  -  j)  +  (U(x  +  |)  -  IJ(x  -  -)]  +  h(-x  -  |)  U(-x  - 

(50) 

where  h(x)  is  an  unknown  function  representing  the  potential  distribution 
at  y  =  h  +  s  for  |xj  >  w/2.  Also,  U(x)  is  the  unit  step  function  with 
the  usual  definition  u(x)  =  1,  x  >  1  and  zero,  otherwise.  The  Fourier 
transform  of  4>(x,  h+s)  can  now  be  written  as 

?(B,  h+s)  =  2  -S-— Bsw/2)  +  eJSw/2  H+(S)  +  e~J6w/2H+(-$) 

(51) 

where  H+(B)  is  the  Fourier  transform  .of  h(x)U(x).  The  subscript  + 
has  been  introduced  to  indicate  that  H+(B)  is  regular  ir.  the  upper  half 
of  the  complex  6-plane.  We  also  note  that  in  view  of  the  symmetry 
in  the  space  domain  <J>  is  symmetric  in  6  in  the  transform  domain. 
Substituting  (51)  in  (25)  we  get 
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(52) 


SW.PJBI  =  2  sin^W21+ej6w/2  +  e~jSw/2  R  (_g) 

o 

which  is  the  desired  equation  to  be  solved  for  the  two  unknowns  p(6) 
and  H+(8).  This  can  be  done  by  an  extension  of  the  Wiener-Hopf 
technique  introduced  by  Mictra  and  Lee  (1971).  In  order  to  extract 
the  solution  for  the  two  unknowns  from  a  single  equation,  we  must  make 
use  of  certain  additional  information  available  to  us.  First,  p(6), 
which  is  the  Fourier  transform  of  a  function  with  a  finite  support 
in  the  x  domain,must  be  regular  in  the  entire  finite  region  of  the 
complex  6-plane.  Second,  the  other  unknown  H+(6)  must  have  all  of  its 
singularities  only  in  the  lower  half  of  the  6-  plane.  We  will  now 
show  how  these  analyticity  properties  of  the  functions  involved  may 
be  put  to  advantageous  use  for  constructing  a  matrix  equation  that 
is  capable  of  yielding  accurate  solutions  in  an  efficient  manner. 

As  a  first  step,  g(6)  is  rewritten  as  follows 


where  K  is  a  known  constant;  I  (o  )  0,  I  (6  )  *•  0,  I  (v  )  '  0,  :  a 

o  mnmnnn  n 

\ 

are  the  pole^;  and  :yn  are  zeros  of  g(B)e  Since  both  p(B)  and 

the  first  term  of  the  right-hand  side  of  (52)  are  regular  in  the  finite  part 
of  the  8-plane,  the  lower  half-plane  singularities  at  6  *  -a  ci  the 
left-hand  side  of  (52)  must  coincide  with  the  poles  of  H+(8)  Likewise, 
tne  singularities  at  6  =  +0^  must  be  coincident  with  the  poles  of  H+(-6). 
Furthermore,  H,(6)  and  H,(-8)  can  have  no  other  singularities. 

T  T 


4L 


Substituting  (53)  in  (52)  and  evaluating  both  sides  of  the 


resultant 

equation  at  B  1  8  , 

II 

B 

m  =  1,  2,  • 

•  •  «  we  obtain 

00 

7 

C'  oo 

n  ,  ,  y 

c 

n 

1  -  > 

m 

(54a) 

L 

n=l 

a  -  6  m 

n  m  n=l 

a  +8 
n  m 

"  ^m 

m  =  1,  2,  •  •  •  » 

00 

r>  CO 

Cn  .  ,  v 

c 

n 

1  ”  ^m 

(54b) 

L 

n=l 

a  -  y  ^m 

n  m  n=l 

a  +  y 
n  m 

jYm 

where 

j8  w 
m 

K  -  e 
m 

'  V 

jY  w 
m 

1  e 

(55) 

The  solutions  of  the  simultaneous  equation  of  the  type  (54)  have 
been  discussed  in  a  number  of  recent  publications  (Mittra  and  Itoh, 
1970,  1972;  Mittra  and  Lee,  1971)  and  will  not  be  elaborated  on. 

The  essential  step  m  deriving  the  solution  entails  the  construction 
of  a  complex  function  f(w)  with  the  following  properties: 


(a)  f(a>)  has  simple  poles  at  a  =  n  =  1,  2,  •  *  •  *>  and  at 


a)  =  0; 


(b)  f (8  )  +  »  r(-S)  -  0 
m  mm 


m  =  1,  2, 


•  •  •  oc 


f(y  )  +  5  f(-Ym)  =  0; 
m  m  m 

(c)  f(u>)  v  | u) |  for  | j  -  >*; 


(d)  the  residue  of  f(ai)  at  w  =  0,  say  R^(0),  is  -  j .  If  such  a 
function  is  available,  then 


c 

n 


=  Rf(«n). 


(56) 
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The  function  f(w)  may  be  expressed  as  (Mittra  and  Itoh,  1970) 


f(w)  =  Kf  (w)P(w) 
o 


f  (w)  =  e^W  f 
0  n=l 


0 

n 


h  -i! 

I  ■r/ 


n  I 


1-2 

o 


(57a) 

i57b) 


where  K  is  a  constant  determinable  by  the  application  of  condition  (d) 
above,  and  P(w)  is 


P(w)  =1+1 


n 


X 

V 


n=l  1  -  —  n=l  1 - 

6  v 

n  n 

F  and  G  are  as  yet  undetermined.  The  factor  eLw  in  (57b),  in  which 
L  is  a  known  constant,  ensures  the  algebraic  nature  of  fQ(w),  viz-, 
fQ(w)  v  1  w |  as  |w|  * 

There  still  remains  the  task  of  determining  the  infinitely  many 

unknown  constants  and  Gn>  However,  the  summation  in  (58)  can  be 

effectively  truncated  at  a  finite  number,  say  n  *  N,  with  N  small, 

because  \  and  i  decrease  exponentially  with  m  as  $  and  L  have  positive 
mm  r  m  m  r 

imaginary  parts.  Hence,  for  n  >  N,  the  condition  (b)  can  be  satisfied 

by  f(B  )  =  f(y  )  =  0.  There  are  now  2N  unknowns  F  and  G  These 
mm  n  n 

may  be  obtained  by  substituting  f (w)  in  the  condition  (b)  followed 

by  the  inversion  of  the  resultant  matrix  equation  for  F  and  G  . 

n  n 

Typically,  highly  accurate  solutions  are  obtained  with  a  matrix  size 
2N  equal  to  10  or  less.  This  is  due  mainly  to  the  fact  that  the 
asymptotic  nature  of  the  function  f (w)  is  incorporated  in  the  solution 
process  making  the  series  representation  a  highly  convergent  one 

Once  the  function  f(w)  has  been  generated,  c^  can  be  obtained  by 
taking  the  residues  of  £(w)  at  w  =  on.  H+(6)  and  p(B)  may  be  readily 


w 


(58) 
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obtained  once  the  coefficients  c  are  known.  The  potential  distribution, 

n 

e.  g.  at  y  =  h,  and  the  charge  distribution  on  the  strip  are  obtained 
by  taking  the  inverse  transforms,  specifically 

»  jo  x 

h(x)  =  -j  l  c  e  n  (59) 

u-1 


The  infinite  series  in  (59),  (60),  or  (61)  converges  very  efficiently 
because  of  the  decaying  exponential  multicative  factors  appearing 
in  the  summation  terms.  It  is  interesting  to  note  that  the  first 
two  terms  in  (61)  correspond  to  the  parallel  plate  capacitance  while 
the  infinite  summation  represents  the  edge  capacitance.  Thus,  the 
fringing  effect  is  accurately  estimated  by  the  application  of  this 
method . 
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Furthermore,  a  knowledge  of  the  rate  of  fall-off  of  the  potential 
away  from  the  center  strip  is  useful  in  estimating  the  coupling 
between  two  adjacent  lines  on  the  same  substrate. 


/, 

III.  WAVE  THEORY  ANALYSIS  OF  MICROSTRIP  LINES 

3.1  Preliminary  Discussion 

In  the  last  section  a  number  of  quasi-TEM  techniques  were  presented 
for  calculating  the  characteristic  parameters  of  ruicrostrip  lines. 

However,  strictly  speaking,  since  the  wave  velocity  in  the  dielectric 
is  different  from  that  in  free  space,  it  is  not  possible  to  support 
a  TEM  mode  in  the  structure.  In  fact,  one  can  show  that  not  even 
pure  TE  or  TM  modes  can  exist  in  the  structure  but  that  it  can 
only  support  a  hybrid  mode  (i.  e. ,  one  in  which  both  the  longitudinal 
electric  and  magnetic  components  are  nonzero)  This  can  be  seen 

rather  easily  from  the  following  consideration  Consider  the  cross 
section  of  a  shielded  microstrip  line  shown  in  Figure  20.  If  the 
center  strip  is  removed  from  the  waveguide,  it  reduces  to  a  partially 
filled  guide  that  can  support  longitudinal  section  electric  (LEE)  or 
longitudinal  section  magnetic  (LSM)  types  of  mode  (see  for  instance 
Collin,  1960),  but  not  a  pure  TE  or  TM  mode.  The  insertion  of  the 
center  strip  in  the  waveguide  causes  currents  to  flow  in  both  the  x- 
and  z-  directions  on  this  strip.  These,  in  turn,  serve  to  couple  the 
LSE  and  LSM  modes  so  that  the  final  mode  configuration  in  the  shielded 
microstrip  line  is  hybrid  in  nature. 

Although  this  fact  has  been  ?  ecognized  for  quite  some  time 
(Deschamps,  1954;  Wu,  1957),  a  rigorous  full-wave  analysis  of  the 
dispersion  characteristics  of  the  microstrip  line  was  not  carried  out 
until  quite  recently.  The  following  paragraphs  will  describe  a 
number  of  these  full-wave  techniques  suitable  for  calculating  the 
dispersion  characteristics  of  microstrip  lines  at  arbitrary  frequencies. 


We  will  begin  with  the  hybrid-mode  analysis  of  the  shielded 
version  of  the  microstrip  line  and  follow  this  with  a  discussion  of 
the  open  microstrip  line.  The  important  steps  in  this  analysis  are 

(i)  Representation  of  the  field  components  in  terms  of  E- 
and  H-  type  of  scalar  potentials; 

(ii)  Application  of  boundary  and  interface  conditions  to  the 
field  components; 

(iii)  Derivation  of  a  characteristic  equation  for  the  propagation 
constant  in  the  guide;  and  finally, 

(iv)  Solution  of  the  characteristic  equation  and  computation  of 
the  dispersion  (k  -  3)  diagram. 

3.2  Shielded  Microstrip  Lines 


Figure  20  shows  the  cross  section  of  the  shielded  microstrip 
line.  The  center  strip  is  assumed  to  be  infinitely  thin  and 
perfectly  conducting.  In  addition,  the  dielectric  material  and  the 
metal  shielding  are  assumed  to  be  lossless.  It  is  well  known  that 
the  hybrid  field  components  can  be  expressed  in  terms  of  a  super¬ 
position  of  the  TE  and  TM  fields,  which  are,  in  turn,  derivable  from 
scalar  potentials  and  (Harrington,  1961).  The  total  field 

can  be  written  as  follows: 


hi  ’  j  "S —  Tie><x’  y)  e"j62 


H,i  -  J  ^wc*.  y)  *~ie° 


(62a) 


(62b) 


Eti =  Vi6)  (x>  y)  e"d6z "  T  2  x  Vih)(x’  y>e"j6z<62c> 
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(62d) 


toe . 


Hti  = 


z  x  Vt^e)(X,  y) 


e"j6z  + 


Vt^h)(x,  y)  e'jBz. 


The  parameter  3  in  the  above  equation  is  the  unknown  propagation  constant, 
z  is  the  z-directed  unit  vector, to  is  the  operating  frequency 

ki =  “VwT »  k2 =  *  (63) 

e  is  the  relative  dielectric  constant  of  the  substrate,  and  e  and 
r  o 

are  the  permittivity  and  permeability  of  free  space,  respectively. 

The  subscript  t  in  Equation  (62)  denotes  the  transverse  coordinate 

variables  x  and  y  and  the  superscripts  (e)  and  (h)  are  to  be  associated 

with  the  TM  and  TE  types  of  fields,  respectively.  The  subscripts  i  =  1,  2 

serve  to  designate  the  regions  1  (substrate)  or  2  (air) . 

From  the  symmetry  of  the  structure,  it  is  clear  that  two  orthogonal 

sets  of  modes  exist,  one  of  which  has  a  symmetric  E ^  and  an  antisymmetric 

H  component  (E  even  -  H  odd)  while  the  other  is  characterized  by  E  odd 
z  z  z  z 

H  even.  The  dominant  mode  is  the  lowest  order  E  even  -  H  odd  mode  which 
2  2  2 

approaches  the  quasi-TEM  solution  for  low  frequencies.  In  what  follows 

we  consider  only  the  E  even  -  H  odd  modes,  although  the  methods 

z  z 

presented  here  are  equally  applicable  to  the  other  types  of  modes  as  well. 


3.2.1  Various  methods  for  solving  the  shielded  microstrip  lines 
(a)  Integral  Equation  Method 

The  scalar  potentials  and  satisfy  the  two-dimensional 

wave  equation  in  region  1  and  2  as  well  as  the  requirement  that  the 
total  tangential  electric  field  derived  from  them  vanish  on  the 
metallic  periphery  of  the  box.  An  additional  boundary  condition  is 
that  the  total  tangential  magnetic  fields  vanish  at  the  plane  of 
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symmetry  x  =  0.  In  view  of  the  boundary  conditions  on  the  side 
walls,  it  is  appropriate  to  write 


,(h) 


/(h) 


(e) 

sinh 

r  COS 

A 

k  x 

‘n 

n  3 

n 

,  (e) 

» 

sinh 

a(2) 

(h  - 

A 

y)  cos  k  x 

n 

n 

n 

(h) 

i 

cosh 

(1) 

a  v  sin 

k  x 

n 

n  ' 

n 

s(h) 

n 

cosh 

«<« 

n 

(h  - 

y)  sin  k  x 
1  n 

(64a) 


(64b) 


(64c) 


(64d) 


where  k  =  [n  -  (1/2)]  tt/L, 
n 


t(1)  =  /  k2  +  (32  -  e  k2  , 

n  /  n  r  o 


,<*>  -  /  £2  +  e2  -  r2  , 

n  /  n  o 


with  k  =  w  / e  u  ,the  free-space  wavenumber.  The  coefficients 
o  v  o  o’ 

A^e\  and  B^  are  as  yet  unknown, 

n  n  ’  n  n  3 

The  total  fields  derived  from  and  must  satisfy  the 

interface  conditions  at  y  =  d.  Imposing  these  on  the  symmetry 
conditions  with  respect  to  the  y  axis,  one  arrives  at  the  following 
four  conditions  which  are  mutually  independent. 

(1)  Ezl  "  Ez2  0  <  x  <  L 

(2)  Ex1  =  Ex2  0<x<L 

(3)  (a)  E ^  =  0  0  <  x  <  t 

(b)H=H„  t  <  x  <  L 
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(4) 


0  <  x  <  t 


/ 


/ 


(a) 

Exl 

(b) 

Hzl 

0  <  x  <  L  . 


These  conditions  are  now  imposed  on  the  field  components  derived 
from  (62)  and  (64) .  By  using  continuity  requirements  at  y  =  d  and 
applying  the  conditions  that  E  and  E  vanish  on  the  center  strip, 

Z  X 

one  obtains  a  pair  of  coupled  homogeneous  Fredholm  integral  equations 
of  the  first  kind  (Zysman  and  Varon,  1969).  These  are  given  by 


1  [G^1) (6) 
i=l  n 

CO 

l  [G*j3)  (3) 
_i  n 


f  h,(5)  cos  k  E>dE,  +  G  ($)  /  h0(C)  sin  k  £d£]sin  k  x  =  0 

1  n  n  ^  2  n  n 

o  t 

(65a) 

t  *  (A)  L  * 

I  h,  (5)  cos  k  ?d5  +  G^  '(B)  /  h~(£)  sin  k  £d£)cos  k  x  =  0. 

1  n  n  .  2  n  n 

o  t 

(65b) 


Here  G^i^(3)»i  =  1,  2,  3,  4  are  known  functions  of  $,  and  the 
functions  h^(5)  and  h2(?)  are  unknown  functions  of  ? .  These 
equations  may  be  transformed  into  a  matrix  equation  algorithm.  The  roots 
of  the  determinantal  equation  of  the  resulting  matrix  may  be 
sought  with  the  aid  of  a  digital  computer  to  yield  the  desired  values 
of  6. 


(b)  Singular  Integral  Equation  Method 

The  singular  integral  equation  method  is  known  to  provide 
efficient  solutions  to  a  number  of  diffraction  grating  and  waveguide 
discontinuity  problems,  e.  g.,  an  iris  diaphragm  in  a  waveguide. 
Recently,  it  has  been  shown  that  this  technique  can  be  advantageously 
employed  to  solve  the  problem  of  dispersion  in  a  shielded  microstrip 
line  (Mittra  and  Itoh,  1971).  The  beginning  stages  of  formulation 
in  this  method  are  identical  to  the  one  presented  in  the  previous 
section;  in  fact,  the  starting  point  of  this  method  is  Equation  (64) 


52 


of  the  last  section.  However,  from  that  point  on,  the  singular 
integral  equation  method  deviates  fundamentally  from  the  conventional 
integral  equation  approach.  A  matrix  equation  is  also  derived  in 
this  method  but  it  is  totally  different  in  character  and  much  more 
efficient  for  extracting  numerical  solutions  than  the  one  obtained 
from  (65) . 

The  first  step  in  the  singular  integral  equation  approach  is 
to  employ  the  conditions  (1)  and  (2)  to  express  and  in  terms 

of  and  A^ .  The  imposition  of  conditions  (3)  and  (4)  then 

—(g) 

leads  to  the  following  coupled  equations  for  the  coefficients  A^ 

and  A^ 
n 


l  >  cos  k  x  =  0 
n  n 

n=l 


0  <  x  <  t 


(66a) 


T  A('e)  k  P  (8)  cos  k  x  -  T  AW  k  T  (B)  cos  k  x  =  0 
n  n  n  n  L,  n  n  n  n 

n=l  n=l 


t  <  x  <  L 


(66b) 


y  A^  '  k  sm  k  x  -  T  A'  '  k  sxn  kx  =  0  0  <  x  <  t  (66c) 

n  n  n  L,  n  n  n 

n=l  n=l 


y  A(‘e;  Q  (B)  sin  k  x  -  T  AW  W  (B)  sin  k  x  =  0 
n  xi  n  n  n  n 

n=l  n=l 


(66d) 


where 


t  <  x  L 


A(e)  =  A(e)  sinh  a(1)  d 
n  n  n 


(67a) 


(67b) 


A(h)  -f-J-A(h>  sinha(1)d 
n  g  k  n  n 

n 

and  P  (6),  T  (g),  Q  (g) ,  and  W  (g)  are  expressed  as  follows: 
n  n  n  n 

«(1)  s  -  g2  a(2) 

P  (6)  =  e  —  coth  G^d  +  — - r  —  coth  (h  -  d) 

n  r  k  n  1  -  B2  k 

n  n 


k  1  -  e 

+  6  ~TT\  - 7  coth  '  (h  " 

a(2)l-g2 

n 


(68a) 


k  k 

T  (g)  =  g2  -rr r-  coth  a^d  +  rr  coth  a^2^(h  -  d) 
n  a  (2)  n  ^(2)  n 

-  n  n 


(68b) 


Qn(3) 


wn(g) 


k  1  -  e 

_n_ - _£  coth  a  <■  ;  (h  -  d) 

a  2)  1  -  g2 
n 


— - -Try  coth  a(1)d  +  -ryr  coth  a(2)(h 

1  -  g2  a  1  n  a  2) 

n  n 


(68c) 


-  d)  (68d) 


where  g  =  g/kQ  is  the  normalized  propagation  constant. 

At  this  point  Equation  (66)  may  be  transformed  into  an  infinite 
set  of  homogeneous  simultaneous  equations  for  and  via  the 

conventional  technique  of  taking  a  scalar  product  with  a  complete 
set  of  functions  appropriate  for  the  various  ranges  in  x.  The 
solution  for  g  may  then  be  determined  by  seeking  the  zeros  of  the 
determinant  associated  with  the  above  matrix  equation. 

We  will,  however,  depart  from  this  conventional  procedure  and 
instead  transform  (66)  into  an  auxiliary  set  of  equations  with 
rapid  convergence  properties.  To  this  end,  (66a)  is  first  differentiated 


with  respect  to  x  and  the  resultant  equation  is  substituted  into  (66c) . 
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This  yields 


r  — ( g)  A  A 

l  A'  k  sin  k  x  =  0  0  <  x  <  t 

n  n  n  *■ 


(69a) 


1  A  k  sin  k  x  =  0  0  <  x  <  t  . 

n~l  n 


(69b) 


Similarly,  differentiating  (66d)  with  respect  to  x,  we  have,  after 


some  rearrangement 


r  —(e)  A 

l  An  kn  c°s  k  x  =  f(x)  t  <  x  < 

n=l  n  n  n 


(70a) 


_ A  A 

1  An  kn  cos  knx  a 

n*l  n  n  n 


t  <  x  <  L 


(70b) 


where 


f(x)  =  l  (a  A<e)  +  b  A<h))  cos  k 


mm  mm 


g(x)  =  l  (c  A^e)  +  d  A(h))  cos  k 


mm  mm 


a_  =  k 
m  m 


Pm(6)W(6)  -  t<S)Q^(S)- 
P(B)W(B)  -  T(B)Q(B)  . 


(71a) 


a  T  (B)W(g)  -  T(B)W  (B) 
b  =  k  -2 - - — - - 2 _ 

m  P(B)W(B)  -  T(B)Q(6) 

a  P(g)Qm(8)  -  Pm(B)Q(B) 
c  =  k  - 2 - m _ 

m  P(0)W(B)  -T(8)Q(6) 

a  f  P(g)  w  (B)  -  T  (B)Q(-B) 
d  =  k  1 - 2 - 2 _ 

m  I  P(B)W(B)  -  T(B)Q(B) 


(71b) 


(71c) 


C  7  Id ) 
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The  functions  P(8),  T(8),  Q(8),  and  W(3)  are  the  asymptotic  limits  of  Pm(3), 

T  (8),  Q  (3),  and  W  (3)  as  m  +  Explicitly, 
m  m  m 


e  -  32  _9  1  -  e 

P(3)  =  e  +  — - 37  +  3  - 4 

1-  i  i 


1-3 


1-3 


T  (3)  =  232 


Q(£)  » 


1  -  e 
_ r 

1-62 


£r  -  62 

W(3)  =  — - +  1. 


(72a) 

(72b) 

(72c) 

(72d) 


1  -  8‘ 

It  should  be  pointed  out  that  Equations  (69)  cud  (70)  -ce  similar  to 

those  obtained  in  connection  with  the  quasi-static  formulation  of  the  iris 

discontinuity  problems  in  a  waveguide  (Lewin,  1966) .  Their  most  important 

characteristic  is  that  the  pairs  (69a),  (70a)  and  (69b),  (70b)  are  exactly 

invertible  via  a  singular  integral  equation  approach.  That  is,  it  is 

possible  to  express  the  coefficients  and  in  terms  of  integrals 

involving  functions  in  the  right-hand  side,  viz.  f(x),  g(x).  In  the 

capacitive  discontinuity  problem,  the  right-hand  side  is  known  and  the 

unknown  coefficients  are  determined  in  this  manner.  In  the  present  case 

the  functions  f(x)  and  g(x)  themselves  contain  the  unknowns  A  and 

n 

_/i_\ 

A^  and  the  result  of  the  inversion  of  (69)  and  (70)  is  a  homogeneous 
set  of  equations,  leading  in  turn  to  an  eigenvalue  equation. 

Following  the  standard  technique  of  solving  the  singular  integral 
equations  with  the  requirement  that  the  tangential  E  fields  are  zero  on  the 
strip,  we  obtain  the  following  equations: 


£  (k  <$  -  a  D  -  M  K  )A(e)  -  F  (b  D  +  N  K  )AW  =  0 

u  n  nm  m  nm  m  tv  m  nn  n  n  n 


(h) 


m=l 


P  pm 


m  pm 


m  p  m 


n=l 


n  pn 


n  p  n 


p  =  1,  2, 


(73a) 


l  (-c  D  -  X  K  )A(e)  +  F  (k  6  -  d  D  -  Y  K  )A<'h)  =  0 

m=l  m  m  *  m  n=l  q  qn  n  qn  n  q  n  (73b) 

q  =  1,  2,  •  •  • 

where  is  the  Kronecker  delta,  and  all  other  quantities  are  as  follows: 

L 


pm 


D  =  f 

nm  J 


f  (x)  sin  k  x  dx 
m  n 
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K„  =  r/ 


77X 

2L 


c°s 

-  sin  k  x  dx 
sin  0  n 


r  \  2  nx 

fm(x)  ‘  L  C0S  2L 


m-1 

y  P  sin  q0 
L  q=l  mq 


-  P 


cos  0 


mo  sin  0 


The  relation  between  x  and  0  is  given  by 


cos 


7TX  1.  TTt  .  v  ,  1  »  IT t  ,  ..  - 

—  =  2  (cos  “  ”  D  +  2  ( cos  ~  +  !)  COS  0; 


the  coefficients  P  are 
mq 


m-1 

y  P  cos  q0 
q=0  mq 


cos  k  x 
_ m_ 

A 

cos  k^x 


All  other  quantities  M^,  N  ,  X^,  and  Y^  are  constants  proportional 

to  a  ,  b  ,  c  ,  and  d  .  The  detailed  derivations  of  these  quantities 
m’  m’  m’  m  n 

may  be  found  in  Mittra  and  Itoh  (1971) . 

The  solutions  of  the  determinantal  equation  corresponding  to 

(73)  are  the  desired  values  of  0  —  the  propagation  constant. 

Although  (73)  comprises  a  doubly  infinite  set  of  equations,  we  can 

truncate  the  associated  matrix  to  a  small  size  since  a  ,  b  ,  c  ,  d 

m  nr  m  m 

as  well  as  M  ,  K  ,  X  ,  and  Y  decrease  extremely  rapidly  with  m. 
m’  m  m’  m  1 

This  can  be  seen  by  observing  that  the  asymptotic  behavior  of  these 

coefficients  predicts  an  exponential  decay.  For  instance,  a^  contains 

difference  terms  of  the  type  P(8)  -  P  (6)  and  Q(8)  -  Q  (6)  in  the 

m  m 

numerator.  For  large  m, 


P(6)  -  Pm(6)  -  er 


1  - 


m 


coth  a*^d 
m 


a 


+  e 


1  - 


m 


(2) 

zr—  cothcP\h  -  d)] 
m 

< 

m 


(2m  -  l)rr 
L 


d 


+  exp 


(2m  -  l)n 


(h  -  d) 


{-2er  +  0(nf2)} 


The  behavior  of  Q(8)  -  Q  (8)  is  similar  and  it  follows  that  a 

m  m 

decays  exponentially  for  large  m.  Similar  comment.-,  apply  to  the  other 
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coefficients  as  well. 

As  will  be  demonstrated  shortly,  retaining  only  a  single 
equation  from  each  set  is  sufficient  for  accurate  computation  of 
numerical  results.  Using  the  fact  that  =  0  (Mittra  and  Itoh, 
1971) ,  the  determinantal  equation  of  the  truncated  set  may  be 
explicitly  written  as 

D(6)  =  [kx  -  M1(?)  -  Y1(6)  Kx]  -  N.(?)  X^B)  K*  =  0.  (74) 


Clearly,  this  equation  is  much  easier  to  handle  than  the 
determinantal  equation  of  a  large-order  matrix  that  results  from 
conventional  processing  of  (65) .  Though  the  derivation  of  (74) 
requires  advance  analytical  processing,  this  effort  is  more  than 
compensated  for  by  the  numerical  efficiency  that  results  due  to 
the  simplicity  of  the  characteristic  equation. 

(c)  Fourier  Analysis  Method 

In  this  method,  the  imposition  of  the  interface  conditions  to 
Equation  (64)  is  carried  out  as  follows  (Hornsby  and  Gopinath,  1969a) 
The  continuity  condition  on  E^  for  0  <  x  <  L  is  satisfied  by  choosing 

(k2  -  B2)  A(e)sinh  o(1)d  =  (k2  -  32)  B(e)  sinh  a(2)(h  -  d) . 

1  n  n  l  n  n 

(75) 

The  remaining  interface  conditions  (2)  -  (4)  [see  paragraph  following 


Equation  (64)]  can  be  satisfied  by  equating  the  functions  F^,  F^j 
and  F^  to  zero.  These  functions  are  given  by 


F1(X)  "  <(Hzl  -  Hz2>  D(x)+lfH-D(x)]}  . 

3  y=d 


(76a) 


F2(x)  =  {(Hyl  -  Hy2)  D(x) 


(1  -  D(x)]} 

y=d 


(76b) 
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D(x)]} 

y=d 


(76c) 


F3(x) 


{(Hxi  -  Hx2> 


D(x)  +  E  [ I  - 


where  D(x)  is  a  function  which  is  zero  on  the  strip  and  unity  outside  oi 
the  strip.  The  required  conditions  are  met  by  constructing  a  complete 
set  of  functions  and  requiring  each  F^(x),  i  =1,  2,3  to  be  orthogonal 
to  every  member  of  the  set,  A  suitable  complete  set  of  functions 
may  be  chosen  as 


sin 


mux 


cos 


mux 


L  »  -  L 

and  the  orthogonality  condition  is 


m  =  0,  1,  2, 


/  F  (x) 
-L  1 


sin 


cos 


,  mux  .  ,  A 
(  “  )  dx  =  0 


i  =  1,  2,  3. 


(77) 


In  the  actual  computation  the  infinite  series  in  and  4’^,  given 
by  (64),  are  truncated  at  some  finite  value  of  n,  say  N.  Substituting 
the  expressions  for  F^x)  into  Equation  (77)  and  performing  the 


necessary  differentiations  and  integrations,  there  results  a  homogeneous 
set  of  3  (N  +  1)  / 2  linear  equations  for  the  unknowns  ,  A^1^ ,  and 
B'  .  In  order  for  a  nontrivial  solution  to  exist,  the  determinant 
of  the  coefficient  matrix  must  vanish.  Since  this  matrix  is  a  known 
function  of  8,  the  dispersion  relation  can  be  derived  by  seeking  the 
value  of  8  that  makes  the  determinant  vanish.  In  actual  calculation 
by  Hornsby  and  Gopinath  (1969a),  N  was  chosen  to  be  10  and  20.  The 
deviation  between  these  two  choices  of  N  vas  found  to  be  less  than 
2  per  cent. 

(d)  Finite  Difference  Method 

The  finite-difference  method,  which  was  discussed  in  connection  with 
the  quasi-TEM  solution  of  the  microstrip  line  problem,  can  be  easily 
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extended  to  apply  to  the  shielded  microstrip  line  being  investigated 
in  this  section  (Hornsby  and  Gopinath,  1969b).  As  a  first  step  toward 
applying  the  finite  difference  method  to  this  problem,  the  cross 
section  of  the  microstrip  line  shown  in  Figure  20  is  discretized  into 
a  lattice  form  with  the  separation  between  t.he  adjacent  net— points 
equalling  Ax  and  Ay  in  the  x-  and  y-directions,  respectively.  Next, 
one  defines  the  discretized  potentials  as 


Y(e)  =  4/(e>(mAX  nAy),  4'(h')  =  ¥(h)(mAx,  nAy) 
mn  *  J  *  mn  J 


and  uses  these  definitions  in  the  wave  equations 


^  =  o 

(78a) 

i  -  1,  2 

=  0. 

(78b) 

This  leads  to  the  discretized  equations  of  the  form 

,<e)  .(e) 


A  .V(e) 
x  mn 


2(1  +  R2) 

mn 


\W  _  vj/ 

m+l,n  m-l,n  ra,n+l 


rV6) 

m,n-l 

(79a) 


L?^h)  =  2(1  +  R2)  -  4*^  -  Y'-bJ  -  rVH)  ,  -  rV^  . 

x  mn  mn  m+l,n  m-l,n  m,n+J.  m,n-l 

(79b) 


where 


X ^  =  (k2  -  82)  Ax2  ,  R  =  Ax/Ay. 


The  expressions  in  (79)  are  appropriate  for  the  interior  net-»points. 
The  boundary  conditions  and  the  resulting  equations  must  be  modified 
when  applied  to  the  net-points  that  are  located  on  the  boundary  or 
at  the  interface. 


After  properly  applying  the  discretized  version  of  the  wave  equation 


and  the  boundary  or  interface  conditions  at  all  of  the  net-points,  one 
obtains  a  matrix  equation  of  the  form 

A<j>  =  (80) 

where  d>  is  a  vector  whose  elements  are  Y^  orY^  and  A  is  X.  with 

mn  nm  1 

i  =  1  or  2.  The  coefficient  matrix  A  is  sparse  and  has  the  size 
2(MN  -  1)  x  2(MN  -  1)  where  MAx  =  L,  NAy  =  h.  The  matrix  eigenvalue 
equation  (80)  can  be  solved  by  one  of  several  standard  algorithms, 
e.  g.,  the  relaxation  method. 

As  pointed  out  earlier,  there  are  two  main  advantages  of  this 
method:  (i)  the  algorithm  is  straightforward;  (ii)  it  requires 
a  negligible  amount  of  analytical  preprocessing.  However,  the 
disadvantage  is  that  the  size  of  the  matrix  which  has  to  be  inverted 
is  rather  large.  In  fact,  in  using  this  method,  Hornsby  and  Gopinath 
(1969b)  worked  with  a  matrix  size  of  100  x  100  or  larger,  but  the 
results  they  obtained  were  not  comparable  in  accuracy  to  those  derived 
by  using  a  much  smaller  size  matrix  equation  formulated  by  one  of  the 
other  methods  described  earlier. 

(e)  Mode-Matching  Method 

The  method,  to  be  presented  below,  has  been  developed  by  Kowalski 
and  Pregla  (1971)  and  is  capable  of  handling  the  case  of  the  finitely 
thick  center  strip. 

The  cross  section  of  the  right-hand  half  of  the  microstrip  line  is 
subdivided  into  five  subdivisions,  as  shown  in  Figure  21. 

The  E-  and  H-type  scalar  potentials  Y^  and  Y^  in  the  i^1 
region  can  be  written  as  follows: 
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WALL-^ 


Figure  21.  Cross  section  of  the  right-hand  half  t'i  the  microstrip  line  with 
a  finitely  thick  strip  conductor. 
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1  x  t-\  \ 

l  An  sin  5^*  cos  n*1J(y  -  h) 


(81a) 


,00  _ 


l  Bn  cos  5(1  x  sin  nS(y  -  h) 
n=l  n  n+1  n+1 


(81b) 


N+1 

{Cn  COs5n  (x  ~  L)  cos  *n1)(y  "  h>  +  ^  cos  ^2)  (x  -  L) 


(2) 

♦  cos  (y  -  h)} 


(81c) 


}=1  fDn  Sin  ?n+l  (x  '  L)  sin  nn+l  (y  "  W  +  Bn  sin  5^(x  -  L) 


(2) 

*  sin  iT+1(y  -  h) } 


(Sid) 


«Ti 

Jx  {En  cos  Cn  (x  "  L)  cos  nn2)(y  "  d)  +  Fn  cos  5*2)  (*  “  L) 


•  sin  n^2) (y  -  d)) 


(81e) 


where 


Wn  sin  5n+l(x  -  L)  sin  ^+i<y  -  d>  +  H„  sin  S<2>(*  -  L) 


C0S  nn+l(y  ”  d)) 


=  (k2  -  e2  -  n(1)2)1/2 


(n  -  1)  7t 

h  -  d  -  A 


(81f ) 


(82a) 


(82b) 


.(2)  _  (n  -  l)ir 
’n  ~  L  -  t 

,«)  .  (k2 . 62 .  ?  <2)2,1  n 


(82c) 


(82d) 


ko  is  the  free-space  wavenumber  and  6  is  the  unknown  propagation 
constant  in  the  guide.  The  scalar  potentials  in  regions  4  and  5 
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be  similarly  written  with  coefficients  Q  ,  $  ,  R  ,  R  ,  S  ,  and  T  .  The 

n  n  ti  n  n  n 

wavenumbers  in  the  x-  and  y-directions  are  denoted  by  5^ ,  > 

-(2) 

and  n  »  and  are  defined  as  follows 
n 


-(1)  (n  -  l)7i 
nn  "“I 


(83a) 

(83b) 


^(2)  -  U  k2  -  62  -  5(2)2)1/2. 

n  r  o  n 


(83c) 


fe 


It  should  be  noted  that  in  Equation  (81),  N  of  the  TM  components 
and  N  +  1  of  the  TE  components  have  been  retained.  It  should  also 
be  noted  that  the  tangential  electric  field  derived  from  the  scalar 
potentials  in  region  2  satisfies  the  boundary  conditions  of  the  guide 
periphery  BCD  as  well  as  at  the  point  I  for  any  coefficients  Cn>  •  •  •» 
The  next  step  is  to  match  the  tangential  field  components  across 
the  interfaces  between  the  regions  and  take  inner  products  with  a 
suitable  set  of  orthogonal  basis  functions.  The  appropriate  basis 
functions  are 


sm 


cos 


sm 


cos 


sm 


cos 


(nf’y),  d  +  4  <y<  h,  x  =  t 


(g^x),  y  =  d,  d  +  A,  t  <  x  <  L 


(n^y),  0  <  y  <  d,  x  =  t. 


(84a) 


(84b) 


(84c) 


The  above  procedure  leads  to  a  homogeneous  set  of  equations  for 

the  unknown  coefficients  A  ,  •  •  •  T  .  Since  the  orthogonal  functions 

n  n 

are  identical  to  those  used  in  the  series  expansion  of  the  scalar 
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potentials,  many  of  the  submatrices  in  the  matrix  equation  become 
diagonal.  Hence,  it  is  possible  to  eliminate  certain  sets  of 
coefficients  without  having  to  invert  the  entire  matrix.  After 
eliminating  all  of  the  coefficients  except  C^,  Dn>  Q  ,  and  Rn,  one 
obtains  a  (4N  +  2)  x  (4N  +  2)  matrix  which  is  the  final  equation 
to  be  solved. 

In  the  actual  calculations  by  Kowalski  and  Pregla  (1971),  N  =  10 
was  found  to  give  a  reasonable  accuracy.  The  convergence  of  the  solution 
was  checked  by  increasing  the  matrix  size  up  to  N  =  20. 

3.2.2  Numerical  results 

In  this  section  we  will  present  a  few  representative  numerical 
results  obtained  by  the  various  techniques  described  in  earlier 
sections.  Refer  to  Figure  22  which  shows  the  numerical  results  obtained 
by  Mittra  and  Itoh  (1971) ,  and  by  Hornsby  and  Gopinath  (1969b) .  The 
effect  of  dispersion  is  evident  at  high  frequencies  where  the  k  -  6 
diagram  deviates  from  the  linear  curve  representing  the  quasi-TEM 
solution.  The  variation  of  the  guide  wavelength  with  respect  to  frequency 
is  plotted  in  Figure  23  for  a  number  of  different  dielectric  materials. 

The  theoretical  results  on  this  curve  have  been  obtained  by  Zysman 
and  Varon  (1969)  and  Mittra  and  Itoh  (1971) .  The  experimental  results 
are  for  an  open  microstrip  line  and  are  included  here  for  convenience 
of  comparison.  It  is  seen  from  Figure  23  that  the  guide  wavelength 
converges  to  the  quasi-TEM  value  at  low  frequencies  and  approaches 
the  wavelength  in  the  substrate  material  as  the  frequencies  increased. 

An  ingenious  empirical  equation  for  the  normalized  wavelength 
of  the  lowest-order  mode  in  the  microstrip  line  has  been  developed 
by  Schneider  (1972).  His  equation  reads 


1 


(85) 


where 


/  e  e  cc 
/  r  eff 


/  £eff  "n  +/  £r 

f2  +  1 
n 


eeff 


e  +  1 


(86a) 


(86b) 


and  Aq  is  the  free-space  wavelength. 

Before  closing  this  section  it  will  be  useful  to  offer  some 
numerical  evidence  that  the  hybrid  mode  solution  approaches  the 
quasi-TEM  limit  as  the  frequency  is  decreased.  Refer  to  Figure  24  which 
shows  the  magnitude  of  the  ratio  of  the  longitudinal  to  the  transverse 
electric  field  components  as  a  function  of  frequency  (Loadholt,  Mittra, 
and  Itoh,  1971).  It  is  evident  that  this  ratio  continues  to  decrease 
as  the  frequencies  decrease  such  that  in  the  low  frequency  .limit 
the  solution  approaches  the  quasi-TEM  form  which  has  zero  longitudinal 
field  components. 


3.3  Open  Microstrip  Line 

Having  discussed  the  shielded  and  boxed  versions  of  the  microstrip 
line, we  will  now  present  a  number  of  techniques  for  analyzing  the 
open  microstrip  line  shown  in  Figure  25.  Even  when  an  open  microstrip 
line  is  placed  in  a  shielded  environment,  the  effect  of  the  enclosure 
may  be  negligible  in  the  event  that  the  walls  are  far  removed  from  the 
center  strip.  In  this  case  the  open  microstrip  line  solution  corresponds 
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to  a  good  approximation  of  the  solution  for  the  enclosed  case.  Of  course, 
the  open  microstrip  line  is  an  important  structure  in  its  own  right;  it  will 
therefore  be  useful  to  discuss  the  details  of  the  solution  for  this  structure 
particularly  since  some  of  the  methods  to  be  described  below  (the  integral 
equation  method,  Galerkin’s  method)  are  sufficiently  different  from  the  ones 
employed  for  the  other  microstrip  structures  considered  earlier. 

3.3.1  Integral  equation  method 


An  integral  equation  formulation  for  this  problem  has  been  given  by 
Denlinger  (1971)  and  the  development  presented  below  will  be  based  on  his 
work.  The  starting  point  in  this  method  is  Equation  (62)  from  which  it  is 
possible  to  represent  all  of  the  field  components  in  terms  of  the  E-  and  H- 
type  scalar  potentials.  However,  in  contrast  to  the  closed  region  problems, 
the  scalar  potentials  and  can  no  longer  be  expressed  in  series 

forms,  as  the  geometry  under  consideration  is  infinite  in  the  x  direction. 
(See  Figure  25.)  It  is  necessary  instead  to  express  the  potential  in  an 
integral  form  as  follows 
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*00 


X 

v: 
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^e)(a,  y) 

e“Jax  da 

(87a) 
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•"Jax  da 
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(a,  y) 

(87b) 
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and  where  i  =  1,  2  designate  the  region  1  (substrate)  and  2  (air), 
respectively.  It  can  easily  be  seen  that  the  transform  potentials  satisfy 
the  following  differential  equations 


dy 

di 

dy' 


2 

Yi 


Yi 


^e)(«,  y)  -  0 


H^h)(a,  y)  =  0 


(88a) 


(88b) 
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where 


2 

Yi  = 

2  2  2 
a  +  B  -  kf 

i  =  1,  2 

(89) 

k1  = 

2 

w  e  v  e  11 

(90a) 

1 

r  r  o  o 

k2  = 

2 

0)  £  JJ  . 

O  0 

(90b) 

The  next  step  is  to  write  the  following  solutions  for  the  transformed 


potential  in  the  various  regions  involved 

?^(a,  y)  =  A(a)  sinh  y^y  (91a) 

y)  =  B(a)  exp[-y2(y  -  d)]  (91b) 

T^(a,  y)  =  C(a)  cosh  y^y  (91c) 

(a»  y)  =  °(o)  exp[-y2(y  -  d)J.  (91d) 


Note  that  the  above  representations  satisfy  the  appropriate  boundary 
conditions  at  y  =  0  and  the  radiation  condition  at  y  ■*■+<*>.  The  unknown 
coefficients  A(a) ,  B(a),  C(a),  and  D (a),  will  be  determined  by  applying 
the  following  continuity  conditions  at  the  interface  y  =  d  in  the 


Fourier  transform  domain. 

Ezl(a,  d)  =  Ez2(a,  d)  (92a) 

Exl(a,  d)  =  Ex2(a,  d)  (92b) 

Hzl(a,  d)  -  Hz2(a,  d)  =  ^(a)  (92c) 

&X1  (« ,  d)  -  d^  =  (.926) 


where  ^  and  are  the  Fourier  transforms  of  unknown  current  components 
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on  the  strip.  Application  of  the  continuity  conditions  (92)  leads 
to  four  simultaneous  equations  from  which  the  coefficients,  A(a) , 
B(a),  C(a),  and  D(a),  can  be  obtained  in  terms  of  the  two  components 
of  current  on  the  strip  and  the  unknown  propagation  constant  3.  The 
expressions  for  the  real  field  quantities  can  then  be  derived  by  the 
use  of  (87)  and  (62). 

We  now  return  to  the  space  domain  (x,  y)  and  apply  the  final 
boundary  condition,  on  the  strip  in  the  form 


Ez2(x,  d)  ■  0 


dF  Hz2(x>  d)  =  0 


[x|  <  w/2. 


(93a) 


(93b) 


The  use  of  these  equations  allows  us  to  represent  and  in  terras 
of  the  coefficients  B(a)  and  D(a)  only. 

oo  2  2 

f  k2  “  e 

Ez2(x,  y)  =  j  I  - g -  B(a)  exp[-Y2(y-d))exp("jax)d<x 


(94a) 

m  2  2 

f  k2  *  6 

Hz2(x,  d)  =  j  /  - -g -  D(a)  expf-v2(y-d)]exp(-jax)da. 

(94b) 


The  coefficients  B(a)  and  D(a)  are,  in  turn,  expressed  in  terms  of 


‘V 

J  (a)  and  J  (a)  as  follows 

X  2 


He)  -  [Fxb22  +  bn)  ixM  +  +  b12\(«) 

*1  '  6 


(95a) 


D<“  -  dat  lFlb21  +  bU>  +  dfe  buV°> 

k.  -  S 

J-  / 


(95b) 
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Our  objective  is  to  derive  a  set  of  equations  for  the  unknowns  J  (a) 

X 

a, 

and  Jz(°0*  This  is  accomplished  by  substituting  (95)  into  (94)  and 
applying  the  boundary  conditions  (93)  to  (94)  resulting  in  the  following 
coupled  integral  equations  for  the  two  current  components. 

00  CO 

I  /  G  (a,  0)  I  (a)  e-jaxda  +  I  f  G19(a,  6)  I  (a)  e"jaxda  =  0 
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where 


-w/2  <  x  <  w/2 


(96b) 
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x  xo  x 


J  (a)  =  I  I  (a), 
z  zo  z 


(98a) 
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For  simplicity,  we  have  used  a  one-term  approximation  for  the  unknown 
current  distributions  and  have  written  them  in  the  form  of  (98a)  and 


(98b).  We  have  also  assumed  that  I  and  I  have  known  forms, and 

x  z  * 

the  only  unknowns  in  their  representation  are  the  amplitude  coefficients 

I  and. I 
xo  zo 

The  unknown  propagation  constant  8  can  now  be  solved  by 
equating  the  determinant  of  the  coefficient  matrix  for  the  unknowns 
I  and  I  to  zero.  Obviously,  within  the  one-term  approximation 
being  used  here  the  results  are  critically  dependent  on  the  choice  of 
the  forms  of  Che  current  distribution.  The  distribution  chosen  by 
Denlinger  (1971)  is  in  the  following  form. 


^  ,  s  l  (  2h  3f (aw)2  -  81  ,  aw  .  [(aw)2  -  12]  (aw) 

I  (a)  -  — / - x  +  o  cos  (  ~xr  )  +  x  sin  ^ 

2  m\(c«)3  (a«)3  2  (0„)2  2 


(99a) 


^  fsin  G^a)  sin  G2(a)  cos  0.4  G3(a)  -  cos  0.5  G3(a) 

Ix(a)  =  G^O  G2(a)  +  ^ 


G2(a) 


G3(a) 


cos  0.4  G,(a)  -  cos  0.5  G,(a] 

+ - 5^> - — i99b) 


where 


G1(a),  G2(a)  =  0.4  (  q-j  ±  aw) 


G3(a),  G^(a)  =  aw  ±  5n. 


Equations  (99)  are  the  Fourier  transforms  of 
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otherwise 


x  [  <0.8  w/2 


Ix(x)  =  C0S  07 


0.8  w/2  <  x  <  w/2 


otherwise 


'  / 


It  may  be  noted  that  the  trial  distribution  for  vhe  current  I2(x)  in 

Denlinger's  work  is  identical  to  that  used  by  Yamashita  (1968)  for 

calculating  the  quasi-TEM  line  capacitance. 

When  the  strip  width  is  small  as  compared  to  the  wavelength,  i.  e., 

w/A  <  0.1,  and  when  we  are  considering  the  lowest  order  hybrid  mode,  it 

is  a  good  approximation  to  neglect  the  transverse  current  and  to 

satisfy  (93)  only  at  the  center  of  the  strip.  The  above  implies  that 

the  boundary  condition  on  the  strip  is  imposed  only  at  the  midpoint 

x  =  0  instead  of  the  entire  range  -w/2  <  x  <  w/2.  This  artifice  allows 

one  to  reduce  the  coupled  pair  of  integral  equations  (96)  into  a 

single  integral  equation  containing  the  unknown  propagation  constant 

6  as  a  parameter.  Note  that  the  above  approximation  also  implies  the 

setting  of  I  equal  to  z.ero  and  x  =  0  in  (96a).  In  addition,  Equation 
xo 

*\/ 

(96H)  now  becomes  trivial  by  virtue  of  the  fact  that  G01I  is  an  odd 

*Lju  Z 

function  of  a. 

This  approximation  is  made  because  a  complete  solution  of  the 
coupled  pair  of  integral  equations  (96)  requires  a  rather  large 
amount  of  computer  time.  Nevertheless  some  test  calculations  have 
been  carried  out  using  the  complete  equations  in  (96),  and  it  has 
been  found  that  the  approximate  solution  obtained  by  letting  I  equal 
zero  is  in  excellent  agreement  with  the  more  exact  solution,  in  the 
zero  frequency  limit  as  well  as  in  the  finite  but  moderate  frequency 
range.  This  is  true  as  long  as  the  normalized  strip  with  w/A  is  less 
than  0.1. 

Fl?.:re  26  shows  the  dispersion  diagram  computed  by  the  single 

2 

integral  equation  approximation.  The  definition  of  is  (Aq/A  ) 
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FREQUENCY  (SKi) 


Figure  26.  Effective  dielectric  constant  versus  frequency 
w/d  =  0.96;  d  =  0.317  cm. 


behavior  agrees  quite  well  with  the  experimental  results  of  Hartwig, 

Massd  and  Pucel  (1968)  even  beyond  the  frequency  range  over  which  the 
transverse  current  is  negligible.  The  curves  also  show  that  the 
effective  dielectric  constant  approaches  in  the  high  frequency 
limit,  indicating  that  at  high  frequencies  all  of  the  energy  tends  to 
be  confined  in  the  dielectric  substrate.  Recall  that  this  phenomenon 
is  similar  to  the  one  observed  in  the  case  of  the  shielded  microstrip 
line  discussed  in  the  previous  section. 

3.3.2  Galerkin's  method  in  the  Fourier  transform  domain 

In  the  previous  section  we  discussed  a  solution  of  the  coupled 
integral  equations  for  the  open  microstrip  line.  An  inherent  step 
in  the  solution  was  the  assumption  of  a  suitable  form  for  the  components  of  the 
current  on  the  center  strip.  Obviously,  the  result  obtained  by  this 
method  is  critically  dependent  upon  the  accuracy  of  the  assumed  forms 
of  the  distribution  for  the  currents.  A  method  is  now  presented  for 
circumventing  this  difficulty  and  systematically  solving  for  the  currents 
to  the  desired  degree  of  accuracy.  The  method  is  basically  a  modification 
of  Galerkin’s  appraoch  adapted  for  application  in  the  Fourier  transform 
domain,  developed  by  Itoh  and  Mittra  (1971a).  It  is  quite  similar  to 
the  one  discussed  earlier  in  Section  2.4  xn  connection  with  th<s  derivation 
of  line  capacitance  in  the  quasi-TEM  approximation. 

We  start  with  Equation  (95)  which  relates  the  unknown  coefficients 

f\j 

B(a)  and  D(a)  in  terms  of  the  two  current  components  J  (a)  and  J  (a). 

Rather  than  applying  the  boundary  condition  (93)  on  the  strip  in  the 
space  domain  as  was  done  in  the  previous  section,  we  impose  this  cc ndition 
in  the  Fourier  transform  domain  instead.  As  a  first  step  we  rewrite 
(93)  as 
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-w/2  <  x  <  w/2 


2  2 
-  6 

j  - g - u(x),  otherwise 


(100a) 


57  h22(x'  d)  * 


-w/2  <  x  <  w/2 


2  2 
k2  "  6 

i  - - - v(x) , otherwise 

'  p 


(100b) 


The  Fourier  transforms  of  the  spatial  functions  in  (100)  are 

2  2 
k  -  $ 

Ez2(a,  d)  =  j  - -  [^(a)  +  fr2(a)] 

k2  -  B2 

^  Hz2(a,  d)  =  j  — g -  t^(a)  +  V2(a)] 
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Using  the  expressions  given  by  (62)  and  (95)  we  obtain  the  following 
coupled  equations  for  the  two  current  components 


Gll(a’  6)^x(a)  +  G12(a>  3)Va)  =  Ux(a)  +  U2(«)  (102a) 

G21(a’  6)^(a)  +  G22(a*  6^z<a)  =  Va)  +  Va)  (102b) 
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if 
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where  G^,  g21>  and  G22  are  8iven  by  (97).  Note  that  the  two 

equations  in  (102)  actually  contain  six  unknowns  altogether.  However, 

four  of  the  unknowns,  viz.,  t^,  U2>  V^,  and  ,  can  be  eliminated  from 

these  equations  by  using  a  technique  that  was  outlined  in  Section  2.4. 

As  a  first  step  we  expand  J  and  J  in  terms  of  known  basis 

x  z 

functions  J  and  J  as  follows: 
xn  zn 


n  M 

J  (a)  =  l 
x 

c  J  (a) 
n  xn 

(103a) 

n=l 

N 

3fz(a)  -  l 

"  n=l 

dn  («>. 

n  zn 

(103b) 

f\j  f\j 

The  basis  functions  J  (a)  and  J  (a)  must  be  chosen  such  that  their 

xn  zn 

inverse  Fourier  transforms  are  nonzero  only  on  the  strip  -w/2  <  x  <  w/2. 
After  substituting  (103)  into  (102)  we  take  the  inner  products  with 
the  basis  functions  Jxn  and  J  for  different  values  of  n.  This  yields 
the  matrix  equation 

0  m  —  1,  2,  •  •  •  M 

(104a) 

0  m  -  1,  2,  •  •  •  N. 

(104b) 

^xn(»)  d«  (105a) 

^zn(a)  da  (105b) 
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An  application  of  Parseval's  theorem  will  verify  that  the  right-hand 
sides  of  (102)  are  indeed  eliminated  by  this  procedure.  Using  this 
theorem  we  can  show  that 


/ 

—CO 


J  (a)  [fiL  (a)  +  lL(a)]  da 
xm  i  Z 


J  (x) 

xm 
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3  0«2  " 


a2) 


Ez2(x* 


d) 


dx  =  0 


The  above  relation  is  true  since  J  (x) ,  the  inverse  transform  of  ^  (a), 

xm  xm  ’ 

and  E  „(x,  d)  are  nonzero  in  the  complementary  regions  of  x. 
z  L 

The  next  step  is  to  solve  the  simultaneous  equation  (104)  for 
the  propagation  constant  6,  by  setting  the  determinant  of  this  set 
of  equations  equal  to  zero  and  seeking  the  roots  of  the  resulting 
equation. 

It  may  be  of  interest  to  note  that  the  technique  just  described 
is  useful  for  solving  the  problem  of  the  slot  line  (Itoh  and  Mittra, 
1971b)  which  is  another  useful  transmission  line  structure  used  in 
microwave  integrated  circuits. 
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IV.  HIGHER-ORDER  MODES 


'/ 


Until  now  we  have  restricted  ourselves  to  the  discussion  of  the 
dominant  mode  in  the  microstrip  transmission  line.  In  this  section 
we  will  briefly  take  up  the  subject  of  higher-order  modes  that  can 
propagate  in  these  lines.  Although  much  has  been  written  on  the 
characteristics  of  the  dominant  mode,  relatively  little  has  been 
reported  on  the  subject  of  hig’ner-order  modes. 

Let  us  first  consider  the  open  microstrip  line.  Since  the 
geometry  of  this  structure  is  an  open  one,  its  modal  spectrum  is 
comprised  of  two  parts  —  discrete  and  continuous.  The  discrete  'part 
of  the  spectrum  corresponds  to  a  finite  number  of  propagating  modes 
which  are  essentially  surface-wave-type  in  nature.  Since  the 
propagation  constant  g  for  these  modes  is  real  and  greater  than  the 
free-space  wavenumber  kQ,  they  are  referred  to  as  slow  waves.  In 
contrast,  the  continuous  spectrum  is  associated  with  the  radiated 
field  or  the  fast  waves.  As  the  frequency  of  operation  is  increased, 
leaky  modes  or  pseudomodes  which  are  alternate  representations  of 
portions  of  the  continuous  spectrum,  can  exist  simultaneously  with  the 
surface-wave-type  modes.  These  leaky  modes  represent  radiation  losses, 

tf 

and  the  use  of  the  microstrip  line  becomes  restricted  when  these 
pseudomodes  appear. 

In  addition  to  the  leaky  modes,  there  are  higher-order,  surface-wave- 
type  modes  that  also  appear  as  the  frequency  of  operation  is  increased. 
The  occurrence  of  these  higher-order,  surface-wave-type  modes  may 
be  heuristically  explained  as  follows.  If  the  center  strip  of  the 
open  microstrip  lines  structure  is  removed  the  result  is  a  conventional, 
dielectric  slab,  surface-wave  line  supported  by  a  conducting  sheet. 
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Such  a  surface-wave  line  supports  the  even  TM  and  odd  TE  type  of 
mode  (Collin,  1960).  The  insertion  of  the  center  strip  perturbs 
these  surface  waves,  but  their  modified  versions  remain  similar  in 
character  to  the  dielectric  slab  modes.  The  existence  of  these 
higher-order  modes  is  undesirable  and  the  use  of  the  microstrip  line 
is  typically  restricted  to  frequencies  below  f^,  the  cutoff 
frequency  of  the  TE^  surface  wave  (Denlinger,  1971), which  is  given 
by  the  formula 

f  = - 2HT>  c  =  velocity  of  light  in  vacuum. 

C  4d  Jt~-  1 

Let  us  now  turn  to  the  problem  of  higher-order  modes  in  a  shielded 
microstrip  line  (Figure  20).  As  expected,  this  closed  waveguide 
structure  only  supports  an  infinite  number  of  discrete  modes  and 
the  propagation  constants  for  these  modes  are  obtained  by  solving  the 
characteristic  equation  derived  in  Section  3.2. 

Figure  27  shows  typical  plots  of  the  value  of  the  determinant  versus 

the  normalized  propagation  constant  6  (Mittra  and  Itob,  1971).  The 

curves  for  increasing  frequencies  are  shown  from  (a)  through  (d) .  It  is 

evident  that  the  number  of  zeros  of  the  determinant  increases  with 

increasing  frequency  indicating  the  appearance  of  higher-order  modes. 

Figure  28  presents  the  dispersion  diagrams  for  both  the  dominant  and 

the  higher-order  modes.  It  should  be  mentioned  that  this  mode  spectrum 

is  not  complete  since  only  the  E  even  -  H  odd  type  of  modes  are 

z  z 

exhibited  in  the  figure.  Calculation  of  the  dispersion  characteristics 
of  higher-order  modes  have  also  been  carried  out  by  Daly  (1971)  and 
Kowalski  arid  Pregla  (1971).  A  comparative  study  of  these  works  reveals 
that  some  disagreement  exists  between  the  results  that  are  reported 
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V. 


LOSSES  IN  MICROSTRIP  LINES 


The  discussion  of  the  microstrip  line  thus  far  has  been  based  on 
the  assumption  that  the  losses  are  negligible.  However,  in  practice, 
losses  are  always  present  and  it  may  be  quite  important  to  obtain  a 
quantitative  estimate  of  these  losses  in  order  to  reliably  design  the 
microstrip  line  circuits.  It  will  therefore  be  pertinent  to  include 
here  brief  descriptions  of  two  methods  that  have  been  found  ujseful  for 
evaluating  these  losses.  We  restrict  our  attention  to  the  dielectric 
and  ohmic  skin  losses  only  and  assume  that  magnetic  losses  in  the 
substrate  are  either  absent  or  negligible. 

The  analysis  will  be  based  on  a  perturbational  approach  which  is 
valid  when  the  loss  per  unit  length  is  small.  Except  for  the  case 
when  the  substrate  material  is  a  semiconductor,  the  assumption  that 
the  losses  in  the  line  are  small  is  certainly  valid  for  most  practical 
microstrip  structures.  The  line  loss  can  be  quantitatively  described 
in  terms  of  an  attenuation  constant  a  defined  as 


P(z)  =  Pq  exp (-2a z) 
P(0)  =  P 

o 


(106) 


where  P(z)  is  the  power  transmitted  along  the  line  at  the  distance 
z  from  the  origin.  Letting  a  =  a^  +  i.  e.,  identifying  the  total 
attenuation  as  the  sum  of  the  contribution  due  to  dielectric  and 
ohmic  losses,  we  obtain  the  following  relationships 


a 


dP/dz  Pc  *  Pd 
2P(z)  %  2P(z) 


(nepers/m) 


(107a) 


and 
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The  perturbational  analysis  for  loss  calculation  is  based  upon  the 
assumption  that  the  field  distribution  in  the  structure  is  not 
altered  due  to  the  presence  of  losses.  Assuming  further  that  the 
dominant  mode  is  quasi-TEM,  the  following  formula  may  be  written  for 
the  two  components  of  a(Yamashita  and  Atsuki,  1970) 

If  a ,  (Vij>)2dxdy 

a  =  - - - 5 -  (108a) 

2  //  ve(V4>)  dxdy 

/  R  i  ZdZ 

a  =  - — - 5 -  (108b) 

c  2  //  ve(7<J>)  dxdy 

i  =  p  v  (108c) 

s  s 

where  the  double  integrals  are  defined  over  the  cross  section  and  the 
line  integral  is  taken  around  the  center  strip  and  along  the  ground 
conductor  surface.  The  various  quantities  appearing  in  (108)  are 
pg  =  charge  distribution  on  the  conductors 
4>  =  potential  distributions 

v  =  phase  velocity  (equal  to  the  group  velocity) 

=  conductivity  of  dielectric  material 

R  =  v4)p/2o„  =  surface  resistance  of  conductors, 
s 

The  quantities  <j»,  pg,  and  v  in  the  above  expression  are  calculated  by 
one  of  the  methods  presented  in  Section  II  with  the  assumption  that 
there  is  no  loss  in  the  microstrip  line. 
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Also,  it  has  been  assumed  in  writing  (108)  that  and  Rg  are 
sufficiently  small  (od  «  we;  and  Rg  <<  Z,  the  characteristic  impedance). 

An  alternative  approach  has  been  followed  by  Pucel,  Masse,  and 
Hartwig  (1968)  who  have  developed  formulas  that  express  ac  and  in 
terms  of  structural  parameters  of  the  guide  and  the  filling  factor 
introduced  by  Wheeler  (1965).  The  expression  for  a^,  the  contribution 
due  to  dielectric  losses,  is  given  by 

a ,  ^  4.34  — 9 — •  o  (dB/cm)  (109) 

"  v  /e  CS  o 

/  eff 

where  q  is  the  filling  factor  and  is  the  effective  dielectric 

constant.  These  quantities  were  defined  earlier  in  Section  2.2.  The 
factor  4.34  represents  the  conversion  of  nepers  into  decibels. 

Equation  (109)  is  applicable  to  an  open  microstrip  line  with  a 
single  layer  dielectric  substrate  below  the  center  strip  and  free 
space  above. 

To  obtain  the  ohmic  attenuation  constant  a  ,  Pucel  et  al.  used 

c 

a  technique  based  on  the  so-called  "incremental  inductance  rule" 

(Wheeler,  1942).  This  rule  expresses  the  series  surface  resistance 

R  per  unit  length  in  terms  of  that  part  of  the  total  inductance 
s 

per  unit  length  which  is  attributable  to  the  skin  effect,  i.  e.,  the 
inductance  produced  by  the  magnetic  field  within  the  conductors. 

It  is  well  known  that  the  surface  impedance 

Zs  =  Rs  +  jXs  (110) 

has  a  real  part  Rg  (surface  resistance)  which  is  equal  to  the 

imaging  part  X  ,  where 
s 
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X  =  U>L. . 
s  1 


(111) 


According  to  Wheeler,  can  be  inferred  from  the  external  inductance 
L  per  unit  length  as  the  incremental  increase  in  L  caused  by  ,-tu 
incremental  recession  of  all  metallic  walls  carrying  a  skin  current 
(see  Figure  29).  The  amount  of  recession  is  equal  to  half  the  skin 
depth  ,5  =  /2/u)uo  .  An  assumption  underlying  this  rule  Is  that  the 

radius  of  curvature  and  the  thickness  of  the  conductors  exposed  to 
the  electric  field  be  greater  than  the  skin  depth  —  pre; '-.rafcly  several 
skin  depths.  According  to  Wheeler,  we  have 


-  r  ^jn  3L 

i  "  1  u  3n  2 

mo  m 


(112a) 


R  -  J-Siik 
s  uo  3nn 


(112b) 


where  the  derivative  8L/8n  denotes  the  derivative  of  L  with  respect 

m 

to  the  incremental  recession  of  wall  m,  n  the  normal  direction 

’  m 

to  this  wall,  and  R  =  wp  6  /2  is  the  surface  resistance 

sm  m  cm 

of  wall  m.  Thus,  from  the  definition  (107c) 


«  =  - T- 

2 | I | 2Z 


o  m 


sm  on 


(113) 


where  Z  is  the  characteristic  impedance  of  the  micrcstrip  line 
calculated  under  the  quasi-TEM  approximation  for  the  lossless  case, 
and  I  is  the  total  current  per  conductor. 

We  assume  that  the  inductance  per  unit  length  for  the  inhomogeneous 
dielectric  case  (microstrip  line)  is  approximately  the  same  as  that 
of  the  unloaded  TEM  line.  This  assumption  implies  that  the  stored 


magnetic  energy  is  not  affected  by  the  presence  of  the  nonmagnetic 
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dielectric  substrate.  This  is  a  reasonable  assumption  as  verified 
experimentally  by  Pucel,  Massd,  and  Hartwig  (1968). 


The  remaining  task  Is  to  derive  the  expression  for  the  inductance 
L  of  the  lossless  line.  From  the  approximate  results  by  Wheeler  (1965) 
we  get 

L  =  ^  +  (  “V  +  •  •  •],  w/h  <  2  (114a) 


where 


L  =  — - - ± _ 

2  ij  +  7in‘2"e€+0’94)1 


w/h  >  2  (114b) 


i  ,  t  „  ,  4nw  . 
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ir  t 


=  w  +  |  £r.  (  ~  +  1) 


w/h  <  l/2ir  (115a) 
(2t/h  <  w/h,  l/2u) 

w/h  >  l/2u .  (115b) 


Assuming  that  the  surface  resistance  of  both  the  center  strip  and 

the  ground  plane  are  identical  to  R  ,  we  obtain  the  following  result 

s 

for  the  ohmic  attenuation  constant 
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where  a  is  in  dB/cm.  From  the  above  expression  it  is  evident 
c 

that  for  a  fixed  characteristic  impedance  which  implies  fixed  t,  w, 
h  ratio  and  e^,  «c  decreases  inversely  with  the  substrate  thickness 
h,  and  increases  with  the  square  root  of  the  frequency.  The  latter 
results  from  the  fact  that  Rg  is  proportional  to  the  square  root  of 
the  frequency. 

Figure  30  shows  some  numerical  results  calculated  from  (116) . 

Figure  31  shows  the  plots  of  the  results  for  a obtained  by  Yamashita 
and  Atuki  (1970)  using  Equation  (108b) .  In  their  paper  a  comparison 
was  made  with  the  curves  in  Figure  30.  Also,  the  dielectric  attenuation 
constant  was  calculated  and  compared  with  the  experimental  results 
by  Hyltin  (1965).  These  are  plotted  in  Figure  32  (Yamashita  and  Atsuki,  1970). 
It  is  evident  that  the  agreement  between  the  theoretical  and  experimental 
results  is  quite  good. 

The  total  attenuation  constant  a  was  measured  by  Pucel,  Masse, 
and  Hartwig  (1968)  for  a  practical  microstrip  line.  The  results 
for  the  case  of  a  rutile  substrate  are  reproduced  in  Figure  33. 

Because  of  the  scatter  in  measured  data  points,  a  curve  was  drawn 
through  the  points  by  "eye-ball"  averaging  where  possible.  The 
agreement  of  the  experimental  data  with  theory  appears  to  be  quite 
good.  The  sharp  upturn  in  Figure  33c  between  5  and  6  GHz  was  claimed 
to  be  caused  by  the  excitation  of  the  TE^  surface  wave  which  produces  ’ 
a  loss  through  propagation  of  energy  out  of  the  edges  of  tne  substrate. 

This  mode  has  a  cutoff  frequency  £  =  c/(4h)»''er  -  1  v  5.8  GHz. 

Indeed,  Pucel  et  al.  observed  radiation  from  the  sides  of  the  slab  in 
and  about  this  frequency  range. 
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Figure  30.  Theoretical  conductor  attenuation  factor  of  microstrip  as  a 
function  of  w/h,  calculated  by  Pucel  et  al.  (1968). 
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Figure  31.  The  attenuation  constant  due  to  the  conductor  surface  resistance 
of  the  microstrip  line.  [The  other  theory  results  are  from 
Pucel  et  al.  (1968).] 


Figure  32.  The  attenuation  constant  versus  the  conductivity  of  the 
dielectric  substrate  of  the  microstrip  line,  e*  =  11.7. 
[The  experimental  results  are  those  of  Hyltin  (1965).] 
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VI.  RADIATION  AND  END  LOADING 


In  the  previou  sections  we  have  concentrated  our  attention  on 
infinitely  long  microstrip  lines.  In  practice,  however,  the  line 
sections  are  necessarily  finite  in  extent.  Thus,  it  is  necessary 
in  practical  designs  to  evaluate  the  effects  of  discontinuities  or 
junctions  introduced  in  uniform  microstrip  lines.  Typical  discontinuities 
of  interest  are  T-junctions,  gap  in  the  center  of  the  microstrip, 
truncated  microstrip  sections,  etc.  In  this  section  we  discuss  the 
open-ended  microstrip  stub  as  a  typical  example  of  such  junction 
problems. 

The  open-ended  microstrip  stub  is  frequently  used  as  a  component 
for  a  filter  or  a  matching  network.  An  ideal  lossless  stub,  with  a 
true  open  circuit  at  the  end, appears  as  a  pure  susceptance  at  the 
input  junction  where  it  is  connected  to  the  main  line.  However, 
due  to  radid* ion  and  fringing  effects  that  are  always  present,  a  true 
open  cir..  is  never  realized  at  the  end  of  the  stub.  The  effect 
of  the  radiation  can  be  represented  by  a  finite  conductance  G,  whereas 
the  fringing  effects  as  well  as  the  effect  of  higher-order  modes 
generated  at  the  end  of  the  stub  may  be  described  by  a  shunt  susceptance 
B.  The  combined  effect  of  these  is  a  finite,  complex  terminating 
impedance  at  the  end  of  the  stub.  This  and  the  distributed  losses 
in  the  stub  (see  Section  V)  cause  the  input  admittance  as  seen  from 
the  junction  of  the  stub  to  the  main  line  to  be  complex.  In  order  to 
evaluate  this  input  admittance  accurately  in  the  presence  of  the  end 
effects  just  described,  it  is  necessary  to  calculate  the  load 
admittance  at  the  end  of  the  stub.  In  the  following  sections  we  present 
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a  number  of  techniques  for  evaluating  this  load  admittance. 

6.1  Radiation  Conductance 

The  problem  of  radiation  at  the  truncated  end  of  the  microstrip 
line  has  been  considered  by  several  workers  (Lewin,  1960;  Denlinger, 
1969;  Sobol,  1971).  Here  we  follow  a  method  introduced  by  Sobol  (1971) 
for  calculating  the  equivalent  conductance  of  an  open-ended  stub.  The 
geometry  under  consideration  is  shown  in  Figure  34.  In  addition  to 
the  quasi-TEM  approximation  the  following  assumptions  will  be  made: 

1.  h/X  «  1; 

2.  The  field  distortion  at  the  end  of  the  line  is  negligible; 

3.  The  center  strip  is  infinitely  thin; 

4.  Only  the  effect  of  the  dominant  mode  needs  to  be  considered 
at  the  discontinuity. 

It  is  also  assumed  that  the  x  component  of  the  electric  field 
E  is  a  constant,  equal  to  E  in  the  region  |z|  <  w/2  and  |x|  <  h/2 
and  that  it  is  identically  zero  outside  of  this  region.  The  fringing 
effect  of  the  field  is  introduced  via  the  use  of  the  effective 
dielectric  constant  e  ^  defined  by  Wheeler  (1965)  and  discussed 
in  Section  2.2.  The  problem  of  calculating  the  effect  of  radiation 
from  the  aperture  plane  y  =  0  may  be  conveniently  attacked  by  first 
replacing  the  aperture  electric  field  by  a  sheet  of  conceptual  magnetic 
current  bl.  The  magnetic  current  can  in  turn  be  used  to  calculate  the 
magnetic  radiation  vector  L.  The  only  nonzero  component  of  the 
magnetic  current  turns  out  to  be  M^,  and  the  corresponding  magnetic 
radiation  vector  Ln  is  given  by 
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L0  =  -  h  sin  0 


2  cos  6]dz 


(117) 


w/2 

/  Ex  exp[j  f  /sITt 

-w/2 

where  A  is  the  free-space  wavelength  and  0  is  measured  from  the  z- 

axis.  The  total  radiated  power  W,  and  hence  G,  is  given  by 

it  2ir 


W  -7  G(hE  ) 2 
z  o 


sin  0  d<J>  d0 


(118) 


where 


120tt 

reff 


It  should  be  noted  that  we  have  made  the  assumption  that  h/Ao  «  1  in 
deriving  the  above  expression  for  W. 

The  radiation  conductance  G  is  explicitly  written  as 


where 


/• 


„  V  reff  w  / - \ 

G  =  TZT-I(r  /reff 


(119) 


240n 


.  2  ("A reff 

*  sin  “x - 


“  J 


w  cos  0 

/ 


•  3  A 

sin  0 


2 

cos  0 


dO  .(120) 


2 

It  can  be  shown  that  I,  and  hence  G,  varies  as  (w/A  )  for  w/A  «  1 

o  o 

as  w/A  for  w/A  »  1.  For  instance,  for  tt  ft  w/A  much  less  than 
o  o  /  reff  o 

unity 


(121) 


This  approximation  is  quite  accurate  for  (w/A^)  j  ere££  less  than  0.5. 

The  asymptotic  behaviors  of  G  with  resp'rt  to  w/Aq  are  consistent 
with  the  ones  reported  by  Lewin  (1960)  and  Marcuvitz  (1951) . 
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The  driving-point  admittance  of  the  stub  can  be  calculated  in 
terms  of  the  radiation  loading  G,  susceptance  B  due  to  end  effect, 
and  distributed  line  losses  expressed  in  terms  of  Qq (unloaded  Q) . 
The  real  part  g,  of  the  driving-point  admittance,  normalized  to  the 
stub  characteristic  admittance,  is  approximated  by 


81  %  2Q 


2  2  2  2 
6  sec  9  -  tan  6  +  b  tan  9  +  b  (9  sec  9  +  tan  9)  +  g  sec  6 
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’  gd  +  Sr 


(122) 


where  g  and  b  are  the  normalized  values  of  G  and  B,  respectively, 
and  0  is  the  electrical  length  of  the  stub.  The  normalized  conductance 
g^  results  from  the  distributed  losses  and  g^the  radiation  losses. 

The  ratio  of  the  radiation  loss  P  ,/P..  _  ,  for  a  50--S2  line  on 

rad  total 

an  aluminum  substrate  is  plotted  in  Figure  35  as  a  function  of  the 
stub  length.  It  is  evident  that  the  radiation  losses  play  an  important 
role  for  a  short  stub.  The  effective  length  of  the  stub  used  in 
these  calculations  is  obtained  by  extending  it  by  0.4h,  The 
justification  for  this  approximation  appears  in  Section  6.2. 

For  a  quarter-wave  resonator  the  ratio  of  the  radiation  loss 
^rad  anc*  t^ie  distributed  loss  P^  is  given  by 


~  rad 
dist 


(123) 


45rr  1  +  1.6  ~  fz 


ref  fi 


Figure  36  shows  this  ratio  as  a  function  of  frequency  and  c^,  the 
dielectric  constant  of  the  substrate.  The  results  agree  very  well 
with  the  data  presented  by  Denlinger  (1969) . 
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Figure  35.  Ratio  of  radiation  loss  to  total  loss  for  50-0  line  as  a 
function  of  stub  length. 


Figure  36.  Ratio  of  radiation  loss  to  distributed  loss  for  50— £1  line  as 
a  function  of  substrate  dielectric  constant  and  frequency. 
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6.2  Edge  Susceptance 

We  next  turn  to  the  calculation  of  the  lumped  loading  susceptance 
at  the  truncated  edge  of  the  microstrip.  The  susceptance  is  often 
described  in  terms  of  an  edge  capacitance  which  is  represented  by  a 
hypothetical  extension  A£  of  the  microstrip.  The  range  of  A&  is 
approximately  0.2  to  0.5  of  the  substrate  thickness  (Napoli  and  Hughes, 
1971;  James  and  Tse,  1972).  A  simple  theory  for  infinitely  wide  plates 
indicates  an  extension  of  the  length  by  0.44h. 

Recently  a  number  of  attempts  have  been  reported  for  calculating 
the  edge  capacitance  by  solving  the  equivalent  static  problem.  Figure 
37  shows  a  finite  section  of  microstrip  line  of  width  w  and  length  i. 

The  edge  capacitance  for  a  semi-infinite  line  is  calculated  by  first 
obtaining  the  capacitance  for  a  finite  section  of  the  line  and  subtracting 
from  it  the  contribution  of  the  uniform  line.  The  latter  is  equal  to 
the  line  capacitance  for  unit  length  multiplied  by  Z,  the  length  of 
the  line  section.  The  excess  capacitance  is  then  associated  with  the 
fringe  effects  at  the  end.  In  the  following,  two  methods  will  be 
presented.  The  first  method  is  based  on  the  conventional  matrix  equation 
approach.  On  the  other  hand,  in  the  second  method  the  analysis  will  be 
done  in  the  spectral  domain. 

6.2.1  Matrix  method 

We  begin  with  the  three-dimensional  Poisson's  equation 

V2(j>(x,  y,  z)  =  -  —  p (x ,  2 )  6 (y) 

°  (124) 

p(x,  z)  =  0,  | x |  >  w/2,  | z |  >  U2. 

The  Green's  function  G  may  be  defined  as  the  potential  at  (x,  y,  z)  due  to 
the  unit  charge  at  (xo,yQ,zo).  Applying  the  superposition  principle,  it  is 
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€r ,  THICKNESS  h 


Figure  37.  Finite  section  of  microstrip  line. 
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(d  -  b)  +  /(c  -  a)2  +  (d  -  b)2  +  K2 

fn(a,  b;  c,  d)  =  (c  -  a)  l n - -j===. -=====z ::z,n  —  rzzzrr: 

(d  +  Ad  -  b)  +  J (c  -  a)2  +  (d  +  Ad  -  b)2+K^ 

(d  +  Ad  -  b)  +  J  (c+Ac-a) 2+(d+Ad-b) 2+K2 
g^Ca,  b;  c,  d)  =  (c  +  Ac  -  a)  £n - —  —  ■  - - 

(d  -  b)  +  / (c  +  Ac  -  a)2  +  (<?  -  b)2  +  K2 


Cn(a,  b;  c,  d) 


(c  -  a)(d  -  b) _ 

K  /(c  -  a)2  +  (d  -  b)2  +  K2 
n  /  n 


Equation  (126)  is  now  solved  for  the  unknowns  After  letting 

=  1  for  all  i  =  1,  •  •  •  N,  the  total  capacitance  of  the  rectangular 
microstrip  section  of  length  4  is  given  by 

N 

COD  -  E  o..  (129) 

j=l  3 

The  value  of  the  fringe  capacitance  of  the  open-circuited  microstrip 
is  given  by 


C  =  ^  lim  [C(JD  -  £C  ]  (l30) 

ex  2  1  u 

where  is  the  line  capacitance  per  unit  length  of  the  infinitely 

long  microstrip  line,  and  the  factor  1/2  accounts  for  both  ends  of 

the  rectangular  section.  The  limit  appearing  in  Equation  (130)  is 

numerically  computed  as  follows.  The  total  capacitance  for  a  rectangular 

section  is  computed  for  the  sequence  of  increasing  values  of  i,  until 

convergence  is  reached  for  the  computed  value  of  the  excess  capacitance 

C  . 
ex 
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easily  shown  that  /2 

f 


4>(x,  y,  z)  =  /  j  G(x,  y,  z;  xq,  0,  zq)  p(xq,  zQ)dxo  dzQ. 

-A/2  -w/2  (125) 

In  the  method  due  to  Farrar  and  Adams  (1971,  1972)  the  charge 

distribution  p(x,  z)  is  calculated  by  numerically  solving 

Equation  (125)  in  the  manner  described  below.  The  finite  strip  is 

subdivided  into  N  subsections  of  elemental  area  AS.,  It  is  assumed 

3 

that  the  charge  distribution  is  uniform  in  a  subsection. 

Equation  (125)  is  then  discretized  under  this  approximation  and 
yields  the  matrix  equation 

jlwv 


N 


(126) 


where  D^_. ,  the  potential  at  subsection  AS^  due  to  a  uniform  charge 

density  of  magnitude  unity  on  AS.,  is  given  by 

3 


D . .  =  G(x.,  0,  z.:  x.,  0,  z.)  Ax.  Az.. 
ij  i  i  3  3  j  3 


(127) 


The  explicit  form  of  has  been  derived  by  using  the  exact  expression 
for  the  potential  due  to  a  uniformly  charged  rectangular  plate,  and 
applying  the  imaging  technique  successively  across  the  dielectric 
boundary  and  the  ground  plate  to  generate  an  infinite  series  of 
images.  The  final  expression  for  D  is 

09  ,n-l,_..n+l 

D  =  y  - - — — - r*  {f  (x.,  y.;  x.,  y.)  +  g  (x.,  y.;  x.,  y.) 

13  2ireo(l  +  er)  n  i’  3’  ■'j  6n  i*  3  * 3 


-1. 


+  f  (y.,  x.;  y.,  x.)  +  g  (y.,  x . ;  y.,  x.)  -  K  [tan  £  (x. ,y . ;x. ,y .) 
n  Jx  1  '3  3  °n i  3  rr  n  i,Ji’  3  'j 


+  tan’*1Cn^xi>  ^i5  xj  +  Axj»  yj  +  Ay.j>  "  tan-15n(x.,  yi;  x  ,  y^+Ay^) 


-  tan"1Cn(xi,  yi;  x^  +  Ax  ,  y  )]} 


(128) 


where 
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6.2.2  Fourier  transform  method 


An  alternative  approach  to,  the  edge  capacitance  may  be  developed 
by  extending  the  Fourier  transform  technique  presented  in  Section  2,4 
(Itoh,  .Mittra,  and  Ward,  1972).  The  extension  is  necessary  since  we 
are  now  dealing  with  a  two-dimensional  problem  as  opposed  to  a  one¬ 
dimensional  case  discussed  in  Section  2.4, 

The  method  proceeds  by  defining  the  two-dimensional  Fourier 
transform  of  the  potential  function  via  the  equation 

OS  03 


'V. 

<K«» 


y,  6) -  =  J*  j  <K»>  y,  zj  ed djf  c>2>  (131) 


-CO  —CO 


Also,  taking  the  Fourier  transform  of  (124)  we  obtain 


[  •—  -  (a2  +  32)}  <J>(a,  y,  8)  =  -  — -pfr,  8)<5(y)  (132) 

dyz  eo 

where  p  is  the  transform  of  the  charge  distribution 

i/2  w/2 


p(a,  8)  =  f  f  p(x,  z)  e3 («x+8z)  dx  dz. 
-i/2  -w/2 


(133) 


In  view  of  the  boundary  conditions  at  y  =  ~h  and  at  y  =  ",  the  form  of 


the  solution  of  Equation  (132)  is  taken  to  be 

Oi(a,  8)  sinh^/a"  +  82  y  -h  <  y  <  0 

♦  (a*  y»  6)  =  •{ 


(134) 


/  2  2 

B(o,  8)  exp  {-/  *  +8  y)  y  >  0. 


By  proceeding  in  exactly  the  same  manner  as  in  Section  2,4  and  applying 
the  boundary  conditions  on  the  strip,  we  obtain  after  eliminating  A(a,  6) 
and  B(a,  8)  the  following  equation  for  the  transform  of  the  charge 
distribution 

G(a,  8)  p(a,  8)  =  ?(a,  0,  8)  (135) 
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where 


G(a,  B) 


£  ,/ ct^  +  B^  [1  +  e  coth 

o  v  r 


A 


(136) 


2  2 
+  S  h] 


Note  that  G  is  the  transform  of  Green's  function  and  that  the  algebraic 
product  in  the  left-hand  side  of  (135)  corresponds  to  the  surface 
convolution  integral  in  (125).  This  feature  is  very  useful  in  the 
actual  numerical  calculation  because  the  computation  of  the  surface 
convolution  integral  is  a  time  consuming  operation. 

Equation  (135)  is  transformed  into  a  matrix  equation  via  the 
application  of  Galerkin's  method,  which  is  similar  to  that  discussed 
in  Section  2.4.  The  matrix  equation  may  be  written  as 


b  =  l  K 

1Y1  I 


m 


n=l 


mn  n 


m  =  1,  2,  •  •  •  N  (137) 


where 
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mn 


b  = 
m 


CO  00 

/  / 

—CO  —00 

CO  00 

/  / 

—00  —00 


Pm(a,  B)  G(a,  B)  Pn(a,  8)  da  dB 


i/2  w/2 

Pm(a>  B)  ^(a,  0,  B)  da  d8=  jjj-  J  J  p^x, 

-i/2  -w/2 


z)dx  dz. 


The  functions  Pn  are  the  basis  functions  of  p(a,  B)  and  $  is  the 
transform  of  the  potential  on  the  strip,  p^(x,  z)  is  the  inverse 
transform  of  p^fc,  B) .  The  total  charge  on  the  strip  is  given  by 


c(£) 


N/v  /  « •  W  /  £. 

f  r  N 

“  I  CT1  /  /  P_(x>  z)  dx  dz  =  2jt  l  c  b  . 

n=1  n  J  J  n  ,  n  n 

n  Jo/o  _,./o  n=l 


(138) 


Figure  38  shows  the  fringe  capacitance  of  the  open-circuited  micro¬ 
strip  line  calculated  by  the  two  methods  presented  in  this  section. 
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The  values  of  the  fringe  capacitance  can  in  turn  be  used  to  compute 
the  loading  susceptancs  due  to  the  end  discontinuity. 
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.VII .  CONCLUSIONS 


In  the  preceding  sections  we  have  discussed  a  number  of  features 
of  micros trip-type  transmission  lines  for  microwave  integrated 
circuitry.  Approximate  as  well  as  rigorous  analyses  for  the 
characteristics  of  the  infinitely  long  microstrip  lines  have  been 
presented.  Brief  discussion  of  some  practical  problems  encountered 
in  the  actual  design,  such  as  computation  of  losses,  radiation  and 
junction  effects,  etc.,  has  been  included. 

Although  much  has  been  written  on  the  subject  of  microstrip 
lines, there  still  remain  a  number  of  theoretical  as  well  as 
practical  problems  yet  to  be  solved.  Some  of  these  problems  include: 

(i)  comprehensive  analysis  of  higher-order  modes;  (ii)  complex  junction 
problems;  (iil)  radiation  from  a  uniform  section  of  open  microstrip 
line;  and  so  on.  It  is  felt  that  a  simultaneous  development  of  both 
the  theoretical  and  experimental  techniques  will  be  needed  to 
successfully  resolve  these  problems. 
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